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Introduction 


HIRS2, a 20 channel infra-red sounder, and MSU, a 4 channel microwave 
sounder, were first launched on TIRDS^N in November 1978 as an upgraded c^ra- 
ticsial tenperature sounding system. Essentially identical instruments have 
flown on NQAA.-A and NCAA-C and are scheduled bo fly on future operational 
satellites through the eighties. While HIRS2 and MSU were designed primarily 
for the purpose of measuring atmospheric temperature profiles, the observed 
radiances are also sensitive to other meteorological parameters such as sea- 
surface temperature, ground temperature, cloud height and cloud amount, ice 
extent over ocean, sncw cover over land, etc. 

In this report, we describe a physically based research oriented processing 
system for HIRS2/tffiU data, developed at the Goddard Laboratory for Atmospheric 
Sciences (GIAS) , Which provides fields not early of temperature profiles but 
also of all of the above quantities. All parameters are retrieved in a mut- 
ually interacting fashion. Global retrievals have been produced for the period 
Jan 5 - Feb 2, 1979 at a resolution of 125 x 125 km, A preliminary analysis 
of the results shows the accuracy of the atmospheric teonperature retrievals to 
be significantly better than those produced operationally during that period 
using statistical regression techniques. In addition, monthly mean fields of 
retrieved surface tenperature , cloud cover, and ice and snow cover appear to 
be ccnparable to those retrieved from instruments such as AVHRR and SMMR vhich 
were designed primarily to measure these products. 

In all cases, dedicated instruments have to contend with some form of 
"noise", the atmosphere and clouds if one wants to measure surface phenomena, 
the surface and the atmosphere if one wants to measure clouds, etc. While 
the HIRS2/MSU system is not particulary optimized for monitoring any of the 
quantities usually obtained by imagers, a physical retrieval method enabling 
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all quantities to be retrieved simultaneously both minimizes the residual 
effects of atraosEiieric and surface "noise" on the retrieved products and also 
produces fields of all quantities vAiich are ccn^tible with each other. 'Hie 
ability to retrieve climate paraneters from the operational sounding system 
IS all the more significant because this will allow continuous monitoring of 
climate parameters from 1978 on. 

Ihe physical based processing scheme for analysis of the HIRS2/MSU data 
IS an outgrowth of research and experience that began with VTPR data on NOAA-2 
and NOAA-4 [Jastrcw and Halem, 1973; Halem and Sus^ind, 1977] and continued 
with work on HIRS/MSU data from NIMBUS 6 [Halem et al . , 1978] . Earlier ver- 
sions of the processing scheme for the HIRS2/MSU data on TIROS-N are described 
in Susskind and Rosenberg [1979] and Susskind and Rosenfield [1980] . Ihe analy- 
sis procedures are based on the relaxation scheme and cloud filtering methods 
of Chahine [1970, 1974]. 

In the process of performing direct physical solution of the inversion 
problem, atnospheric and surface conditions are found vrfiich meet the require- 
ment that the channel averaged radiances ccnputed from these conditions agree 
with the cbserved radiances to within a specified error limit. Thus, errors 
in corrputation of radiances (hereafter called the forward problon) become a 
serious potential source of noise affecting the accuracy of the method. The 
accuracies obtained by various groups enploying direct physical inversion 
techniques [Jastrcw and Halem, 1973; Halem et al . , 1978; Susskind and Rosen- 
field, 1980] inply the ability to compute expected satellite observations, 
given surface and atmosEheric conditions, with a sufficient accuracy so as not 
to significantly degrade the retrievals. Nevertheless, doubts have been raised 
regarding the current state of the forward problem, that is, the accuracy with 
which our knowledge of atmosE^ieric transmittance functions and our use of the 
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radiative transfer equation assuirptions can reproduce the observed satellite 
radiances for given atmospheric and surface conditions. 

McClatchey [1976] found bias errors of the order of 6°C between radiances 
observed in a number of 15 Vm channels of the DMSP infra-red sounder and those 
calculated using colocated radiosonde observations. In all cases, ccnputed 
radiances were larger than observed ones. Based on his calculations, he w^ 
one of the first to question whether radiative transfer calculations can be 
done in the infra-red region with sufficient accuracy to allow for direct 
physical inversion of infra-red sounder observations. Valovcin (1981) did a 
similar conparision for the 20Vm water vapor sounding channels of the EMSP 
sounder and again found biases in the same sense of the order of 5°C for these 
channels. He suggested further studies of the forward problem to pinpoint 
the sources of these large errors. Weinreb (1979) at NOAA/lJESS showed that 
significant enpirical correctiais had to be made to remove biases between 
observed and coinputed brightness temperatures for the HIRS instrument. No 
indication of the RMS_difference between deserved and confuted radiances was 
given in Vfeinreb (1979) but he iittplies that better accuracy is needed before 
physically based retrievals can ntprove upcai statistically based retrievals. 

In first chapter of this report, we describe the radiative transfer calcu- 
lations performed at GIAS for the HIRS2 and MSU sounders, and show, by direct 
cotparison of satellite observations with those canputed fron colocated 
radiosonde observations, that agreement in both the infra-red and microwave 
channels can be obtained with sufficient accuracy to meet the needs for physi- 
cally based retrievals yielding high accuracy, even in the presence of clouds. 
In the second chapter, a rapid algorithm for ccsiputing the channel averaged 
transmittance functions as a function of tenperature profile, humidity profile, 
ozone profile, and zenith angle of observations is described. Ihis algorithm, 
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which IS used in the analysis of TIRDS^N data, is shown to be sufficiently 
accurate so as not to significantly affect retrieval accuracy. In the third 
chapter, vre describe the theory and application of the processing scheme used 
at GLAS to provide global retrievals (12000 retrievals per day) for the period 
Jan 5 - Feb 2, 1979. In the final chapter, we will show results indicating 
the accuracies of retrieved tenperature profiles for the period as v^ll as 
nonthly mean fields of sea/land surface temperatures, cloud height and cloud 
amount, and microwave surface emissivity at 50.3 GHz, from vrfiich we deduce xce 
cover over ocean and sncw cover ever land. More detailed studies of these 
quantities, as well as a description of inprovements to the processing scheme, 
will be given in future publications. 
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Chapter 1 


An Accurate Radiative Transfer Model for the Direct Calculation of Bri^tness 
Tetperatures for the HIRS2 and MSU Sounding Channels 

1.1 The Radiative Transfer Bquatxon - "The Forweird Problaa’' 

The radiance calculations used in the GLAS pl^ical inversion method are 
based on the form of the radiative transfer equation for ootputing clear column 
radiances R^, for sounding channel x, given by 

Ri = eiBiCTg] Ti(Pg) + (1 - + T^(Pg) 

( 1 . 1 ) 

In P 

+ Pi'Hi-^i'(Ps) + / BiCT(P)] dr dlnP 

In Pg dlriP 

v*iere is the surface emissivity averaged over sounding channel i, B^^CtI 
is the mean Plank blackboc^ function, averaged over channel i, of the tenpera- 
ture T, "riCP) is the nean atmosjiienc transmittance from pressure P to the 
top of the atmosphere and evaluated at 9, the zenith angle of the observation, 
R£+ is an effective atmospheric emission downward flux, [Komfield and Susskind, 
1977], Pj^' IS the directional reflectance of the solar radiation, is the 
channel averaged solair radiance striking the top of the atmosphere, and (Pg) 
is an effective atmospheric transmittance of radiation from the sun, reflected 
frcm the earth's surface to the satellite CSusskind and Rosenfield, 1980]. The 
subscript s refers to the eaxrth's surface. The integral, taken fron the sur- 
face to the satellite pressure P, r^resents the up/elling atmospheric Quitted 
radiation, vhich is a mean value of the bladk-bcx^ function of atmospheric tonr- 
perature weighted by the chcinnel weighting function d^i/dlnP. Table 1 shews 
the channels, centers of the filter functicxis, peaks of the wei^ting functions, 
and other relevant information, for the channels on MSU and. HIRS2. The HIRS2 
and MSU channels used in the GIAS Physical inversion method are designated by *. 
The radiative transfer equation for a plane parallel horaogenecus atmosphere 
expressed by equation (1,1) assumes atmospheric local thermodynamic equilibrium. 
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This assuirption is reasonable for radiance calculations for all channels but 
HIBS2 channel 17/ \riiich contains considerable emission from pressures lower 
than 1 Ttib. For this reason, channel 17 was not included in this study. 


Table 1.1 


HIRS2 and MSU Channels 


Channel 

v(cm“l) 

Pe^ of 

Peak of 





dx/dlnp (mb) 

Bdx/dlnp (mb) 



HI 

668.40 

30 

20 



*H2 

679.20 

60 

50 



H3 

691.10 

100 

100 



*H4 

703.60 

280 

360 



H5 

716.10 

475 

575 

J 


H 6 

732.40 

725 

875 



H7 

748.30 

Suirface 

Surface 



H 8 

897.70 

Window, sensitive 

to water vapor 



H9 

1027.90 

Window, sensitive 

to O 3 



HIO 

1217.10 

lower trcpospheric water vapor 



Hll 

1363.70 

Middle tropospheric water vapor 



H12 

1484.40 

Upper trcpospheric water vapor 



*H13 

2190.40 

Surface 

Surface 



*H14 

2212.60 

650 

Surface 



*H15 

2240.10 

340 

675 



H16 

2276.30 

170 

425 



ffl.7 

2310.70 

15 

2 



*H18 

2512.00 

Window, sensitive 

to solar radiation 



*H19 

2671.80 

Window, sensitive 

to solar radiation 



*M 1 

50.30a 

Window, sensitive 

to surface emissivity 



*M 2 

53.74a 

500 




*M3 

54.96a 

300 




*M4 

57.95a 

70 




a values in CSg 


1.2. Calculation of Atmospheric Transmittances 

The atmospheric transmittance functions x^^CP), contain conponents coming 
from attenuation by discrete lines of absorbing gases, X 3 ^l(P), and also from 
broad-banded ccotinuum absorption features. Ihe ccsiponent coming from discrete 
lines can be calculated by line-by-line calculations according to 

\ 

00 

/ I kj;,(v,z) Cl(z)p(z) dz sec0 

XiL(Pr 6 ) = /dv Fi(v) e - z(P) L (1.2) 
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vrtiere ) is a normalized channel response function, ,z) is the absorp- 
tion coefficient of line L evaluated at the tetrperature and pressure of height 
Z/ cl(z) is the molecular mixing ratio for the gas to vAiich line L belongs, p 
is the density of air, and 0 is the zenith angle of <±)servation. Ihe evalua- 
tion of kL(s) ,z) depends rvDt only on the set of line parameters used [McClatchey 
* et al , , 1973] but also on assutrptions regarding the tenperature dependence of 
the Lorentz half width and the nature of the line shape. Finally, the results 
also depend on the numerical quadrature used in the corputation of equation (1.2). 

All calculations for the HIRS2 channel transmittances were made as m 
Susskind and Searl (1978) using the 1978 version of the AFGL line parameter 
tape [Fothman, 1978] . The HSU transmittance functions were calculated in a 
similar manner, but using a Van-Vleck Vfeiskopf line shape and the overlapping 
line theory given by Rosenkranz (1975) in the case of Og absorption. Conputa- 
tions were done with a 64 level atmosphere and a frequency spacing of .002 cm'“^ 
for the HIRS2 channels and .00006 cm~l for the HSU channels. Filter functions 
for HIRS2 on TIEDS-N were provided by NQAA/NESS. Ihe MSU channels were treated 
as having rectangular response with the specified centers and half-widths. 

The CX >2 line shape was taken to be sub-lorentz as described by Susskind and 
Mo (1978). One significant modification made to the calculations of Susskind 
and Searl (1978) was to include induced emission in the coirputation of the 
tenperature dependence of the line strengths 

S(T) ® “1.439 E"/T [x - e“l * 

I = X _______________ (1-3) 

S(Tg) Q^(T) Qj^ (T) e -1-439 E"/Tg [x - 

where E", Q^. and Qr are defined in McClatchey et al . , (1973). Neglect of the 
induced emission factor, (1 - e“l-439vA)/(l - el-439v/Ts)^ as done in 
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McClatchey et al . , (1973) and Susskind and Searl (1978), decreases the inten- 
sity of lines at low tenperatures relative to high tenperature. For exaitple, 
at V =650 cni“^, the intensity of a line at 220K is underestjjnated relative to 
tive to Its intensity at 300 K by 3%. Such an error has the effect of broad- 
ening the weighting functions of channels sounding the tropopause region. 

Ihe total transmittance function t 3 ^(P) is taken as 

T i(P) = T il(P) T in(P) t ia(P) (1-4) 

where tn/ tvj and xg represent continuum absorption due to N 2 , water vapor, 
and aerosols. Water vapor continuum and Nitrogen continuum absorption are 
treated as in Susskind and Searl (1978). Aerosol absorption and scattering is 
treated according to 

- ka(P)seo9 

Ta(P) = e (1.5) 

where ka(P) is the sum of the aerosol scattering and aerosol absorption opti- 
cal thickness fron pressure P to the top of the atmosphere. The aerosol model 
used IS based on the tropospheric light haze model and calculations of Wang 
and Domoto (1974). Ihe total optical depth of the tropospheric aerosols is 
taken as .01 for nadir viewing in the 15ijm channels. The optical depth falls 
off with a 1.2 km scale height. An additional homogeneous stratospheric aero- 
sol layer between 12 and 20 km was added with total optical depth 30% of that 
of the tropos£iieric layer [Wang, 1975] . Optical depths in the 4.3pm region 
were taken to be a factor of three greater than in the 15pm region, giving a 
total nadir optical depth of .039 for the 4.3pm channels. 

1.3. Need for Rapid Transmittance Algorithms 

Line-by-line transmittance calculations using equation (1.2) take about 
2 minutes CPU (on a 5 MIPS machine) per channel per sounding for the HIRS2 
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channels and is computationally impractical for global analysis of satellite 
data (« 12000 soundings per day) . It is therefore necessary to have a rapid 
algorithm v^ich can reproduce the transmittances/ including their dependence 
on the important variables, the tempera ture-humidity-ozone profile and zenith 
angle of observation, to a desired accuracy. Given an analytical form of the 
approximate trananittance model, line by line calculations are used both to 
generate empirical coefficients and to veri^ the model. An examination of 
equation (1.2) for the special case of monochromatic transmittances 
shows that 

S6C0 

Tv(P,9)=H Tvl(P) (l«6a) 

I 

and Tv(P,e) = Tv(P>,e) Tv(P,P',e) (1.6b) 

vdiere the vertical monochromatic transmittance from pressure P to 

the top of the atmosphere due to absorption from molecular species I, P > P’ , 
and T-y(P', P,0)is the monochromatic transmittance frcsn pressure P to P' at 
zenith angle 6 . Ihese equations are not applicable for channel averaged trans- 
loittances, but they suggest the following form of a rapid algorithm v^ich is 
analogous to equations 6a and 6b; 

_ £ _ _ _ 

Til(Pj^, 6) = n T (P-i ,P-j_j,0) T (Pt,P^^]_,0) t (P^,Pt„j_^®) 

:=1 IF lO iw (1.7) 

where tiL(^£/9)/ the model discrete line transmittance for channel i through 
a path from pressure P^ to the top of the atmos^iiere at zenith angle 6 , is 
given as the product of effective channel averaged transmittances through all 
individual layers from the top of the atmosphere to pressure Pj^. The effec- 
tive layer transmittances, in turn, are modelled as products of three coitpon- 
ents. Hie first, transmittance due to the gases assumed to have 
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fixed concentrations, the second, tiQ/ the effective transmittance due to 
absorption from ozone; and the third, t xW/ effective transmittance 

due to absorption from water vapor. The models for x xF/ lO^ ^ xw 
will be discussed in chapter 2. In all cases, model coefficients are gener- 
ated at a fixed set of angles and the logarithm of the effective layer trans- 
mittance IS taken as being linear in secant 0 . is used in equation (1.4) 

in place of 

1.4, Forward Calculation of Radiances 

In this section, we assume a given tenperature, humidity, and ozone profile 
and compute the transmittance function as described previously. The radiances 
are calculated by numerical integration of equation (1.1) using a 64 level atmo- 
sphere and a Simpson’s quadrature rule. For infra-red channels, the surface 
emissivity ex taken to have the values of .98 over ocean for channels in 
the 12-15pm region and .96 over ocean for channels in the 4.3-3.7ym region [Wolfe, 
1965]. Over land, the emissivity is taken as ,95 and ,85 for the long and short 
wave channels respectively. The effective downward flux Rx+ is modeled according 
to Komfield and Susskind (1977) as 

Ri+ = FiBi(Ts) [l-Tx(Ps)J (1-8) 

The values of Fx found fron this model for the HIRS2 channels are .71, .49, .42, 
and .55 for channels 6, 7, 13, and 14 respectively, and 1.0 for all other channels. 
For the microwave channels, the surface emissivity is calculated from the observa- 
tion of MSU channel 1 according to 

e = Ri - /’D3 t - Rl+xi(Ps) (1.9) 

[Tg - Ri+ ] XI (Pg) 

where R^ is the observed brightness toiperature for MSU channel 1, Rj-l- is the 
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computed downward flux using the teirperature-humidity profile, and '^i(Ps) is 
the surface transmittance for microwave channel 1. The downward flux for micro- 
wave channels is averaged over all angles using an approximation based on mono- 
chromatic transmittances vdiich are linear in the secant of the zenith angle. 

1.5. Accounting for the Effects of Clouds on the Infra-red Observations 

Infra-red radiation is highly attenuated by clouds in the field of view. 

The out-going radiation in cloudy conditions is a complicated function of the 
properties of the clouds as well as the other atmospheric and surface condi- 
tions already discussed. In the analysis of sounding data, we do not attempt 
to conpute the expected outgoing radiation as a function of cloud properties. 
Instead we use a procedure to estimate or "reconstruct" the radiance which 
would have been observed if no clouds were in the field of view by simultaneous 
analysis of observations in two adjacent fields of view [Smith, W. L., 1968; 
Chahine, M, T., 1974], One retrieval is performed in a 125 x 125 Ion area. For 
each area, we separate the many HIRS2 observations (20 x 20 km resolution at 
nadir) in the area into two groups, one set containing the half with the warmest 
observations in the llvtn window, the other set containing the coldest observa- 
tions. Radiances for all spots in each set are averaged together for each chan- 
nels to give representative observations in each "field of view". Ihe observa- 
tions in the two fields of view are used to provide reconstructed clear column 
radiances. The MSU observations assigned to the 125 x 125 km area are those of 
the closest MSU spot to the center of the 125 x 125 km box. No cloud correc- 
tion IS made to the IVBU observations. Atmospheric and surface conditions are 
then found which vhen substituted into equation (1.1) match the reconstructed 
clear column radiances and microwave radiances for the spot. Ife will use the 
same set of observations in the conparison of observed and computed brightness 
temperatures using colocated radiosonde information. 
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The infra-red clear column brightness temperatures are reconstructed by 
making the assunption that the conditicais in both fields of view are otherwise 
identical except for percentage cloud cover, a • Ihe reconstructed clear column 
radiance is then given by 


Ri = + n (Rj.^1 - Ri^2) (1.10) 

vrtiere Rj^^^ is the observation for channel i in the field of view 3 and n is a 
parameter given by a]/(a 2 -ai), where ai and a 2 are the cloud fractions 
in the two fields of view. In analysis of the sounding data, to be described in 
chapter 3, n is determined using observations in HIBS2 channel 13 and MSU 
channel 2 . Ohe micrcwave observations are included in the calculation to mini- 
mize errors in the determination of n due to incorrect estimates of the temp- 
erature profile. For the purposes of this study, v^ere the temperature profile 
is known, the parameter n, vhich is independent of channel, is determined 
according to 


_ TJ 

^ 13 ^13,1 


n = 


R13,1 - R13,2 (1.11) 

where R '^^3 is computed dear-column radiance for channel 13, using the radio- 
sonde temperature profile. 

According to equation (1.10), a value of n = 0 corresponds to field of view 1 
(always taken as the field of view containing the larger radiance in the llym 
window channel) being clear, vhile a value of n = -.5 corresponds to both fields 
of view being clear. 

If the observed channel 13 brightness tenperatures in both fields of view 
differ by less than 1<*C and In I <4, n is set equal to -.5. In this case, both 
fields of view are assumed to be clear and the clear column radiance is given 


12 



by the average of the observations in both fields of view. If a value of n £ 0 
IS ctotained fron equation (1.11) but the difference of the channel 13 brightness 
ten^jeratures is greater than 1®, n is set equal to zero. 

In the GIAS retrieval procedure for the HIRS2/MSU data, atmosjAieric and 
surface prc^rties are found, vdiich, vAien substituted into equation (1.1), match 
the reconstructed clear column radiances, datained from equations (1.10) and 
(1.11). Using the the same procedure, we coipare reconstructed clear column 
radiances with those confuted fron radiosonde conditions colocated with the 
satellite observations. iWo sets of statistics will be presented: cases 

determined to be clear from eq. (1.11), and cases determined to be cloudy. In 
cloudy cases, the computed clear column radiance for channel 13 is artificially 
fixed to agree with the reconstructed clear column radiance by use of equations 
(1.10) and (1.11). In clear conditions, most often represented by n = -.5 as a 
result of homogeneity in both fields of view, a reasonable estimate of the 
accuracy of the calculations for channel 13 is obtained. 

1.6. Conparison of observed and confuted brightness tenperatures 

In this section, to compare the brightness temperatures calculated for 
the HIRS2 and MSU channels from rawinsonde tenperature-humidity profiles with 
observations closely colocated in space and time. The radiosonde profiles 
used in the conparison are a select set of oceanic radiosondes [Phillips ^ 
al. , 1979] vdiich are also colocated in space by 110 km and in tune by 3 hrs to 
the satellite observations. Only those radiosondes reporting tenperatures at 
all mandatory levels from 1000 mb - 30 mb and humidities froti 1000 - 500 mb 
TOre used in the conparison. Hie tenperature above the highest radiosonde 
report was fixed at the values of the stratosjiieric analysis produced by the 
NMC Ipper Air Branch. In cases where the analysis did not extend to 1 mb, the 
temperature profile was extrapolated to climatological 1 mb values. Ccstparisons 


13 



were limited to oceanic regions to minimize the effects of uncertainties in the 
surface ten^rature on the conputed radiances. 

Given the radiosonde tenperature and humidity profile, brightness teitpera- 
tures for the tenperature sounding IR channels on HIRS2 vere calculated accord- 
ing to equation (1.1). 'Ihe teiiperature profile was interpolated to 64 levels assum- 
ing tenperature to be linear in the logarithm of the pressure between the manda- 
tory levels. Ihe surface temperature was taken to be the climatological sea- 
surface tenperature. TSie specific humidity was interpolated between mandatory 
levels assuming a dependence. A specific humidity corresponding to a clima- 
tological water vapor mixing ratio of 3 ppmv was assumed at and above 100 mb. 

Zonally averaged climatological ozaie profiles vrere use to conpute the ozone 
component of the transmittance. 

Two sets of statistics were calculated comparing observed and ccstputed 
radiances. The first sanple represents coirparisons made only in those cases 
where at least one field of view was found to be clear, that is n _< 0. The 
observed radiances in this case correspond to either the average or the warmer 
of the radiances in the two fields of view. Ihe second sample contains only 
partially clou^ cases, n > 0, in vdiich case comparisons were made between 
radiances reconstructed from eqns. (1.10) and (1.11) and canputed clear column 
radiances. The first set of statistics, namely clear cases, is a measure 
primarily of the accuracy of the forward probloon calculation. Tie second set, 
cloudy cases, is also a measure of the nethod’s ability to account for cloud 
effects. 

Ct&ax> Case 

Table 1.2 gives the statistical ccaparisons relating the observed brightness 
temperatures with those computed using colocated radiosonde information for a set 
of 74 clear profiles from Jan-Feb 1979. The first two columns indicate the chan- 
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nel number and characteristic frequency of the channel as shown in Table 1.1. In 
order to minimize the effects of solar radiation, a subset of 21 cases occur ing 
at night were used for comparison in channels 18 and 19, the two shortwave window 
channels and 8, the longwave window channel. While channel 8 is not affected by 
solar radiation, the subset of night cases was taken for consistency of the results 
with those in the two shortwave windows. 

■nie third and fourth cx>lurnns show the untuned mean and standard deviation of 
the differences between ccaiputed and observed brightness teirperatures for each of 
the channels. The average absolute value of the mean difference is .68° C and the 
avarage standard deviation is 1.11° C, The largest standard deviations in the 
infra-red occur in channels 1-3 and 16, vhich are most sensitive to the tenperature 
profile above 20 mb, not well monitored by radiosondes, and in the window channels, 

8, 18, and 19, which are most sensitive to the sea-surface temperature, fixed at 
climatology. Typical deviations of sea-surface temperatures frcM their climatologi- 
cal values can be up to 2° C. The standard deviations in the reiiainder of the chan- 
nels is less than 1°, with the exception of channel 7, which also is affected some- 
what by errors in the climatological sea surface tenperature. The standard devia- 
tions in the microwave channels are slightly larger than in infra-red channels 
sounding coraparable portions of the atmosphere. Perfect "agreement" is forced in 
microwave channel 1 by the surface emissivity determined according to equation (1.9). 
At least part of the randan component of the differences between observed and 
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TABIiE 1.2 


COMPARISC»^ OF OBSERVED AND COMPUTED BRIGEfINESS TEMPERATURES (°C) 
CLEAR CASES JAN - FEB 1979 (74 COLOCATIONS) 

UNTUNED TUNED 




MEAN 

STD EEV 

e 

MEAN 

STD DEV 1 

STD DEV 

STD 1 

CH 

V (cm“l) 

(COMP-OBS) 



(COMP-OBS) 


CBS 

COMP 1 

1 

668.4 

1.55 

2.37 

0 

1.55 

2.37 1 

2.82 

2.83 1 

* 2 

679.2 

-.33 

1.19 

0 

-.33 

1.19 i 

3.86 

3.81 1 

3 

691.1 

1.42 

1.28 

■ 0 

1.42 

1.28 1 

3,60 

2.96 1 

* 4 

703.6 

1.87 

.92 

.12 

.03 

.70 1 

2.32 

2.54 1 

5 

716.1 

.66 

.56 

.03 

.04 

.57 1 

4.31 

4.17 1 

6 

732.4 

.26 

.70 

.015 

-.17 

.70 1 

5.19 

5.03 1 

7 

748.3 

.33 

1.05 

.04 

-.37 

1.05 1 

5.07 

4.74 1 

**8 

897.7 

-.50 

1.86 

0 

o 

LO 

• 

1.86 i 

5.26 

4.25 1 

*13 

2190.4 

1 

• 

to 

o 

.70 

0 

-.20 

.70 1 

6.97 

6,83 1 

*14 

2212.6 

.86 

.68 

.04 

.05 

.65 1 

7.09 

7.19 1 

*15 

2240.1 

,85 

.70 

.035 

in 

o 

• 

.66 1 

6.09 

6.28 ! 

16 

2276.3 

.39 

2.12 

0 

.39 

2.12 1 

2.84 

3.38 1 

**18 

2512.0 

.03 

1.30 

0 

.03 

1.30 i 

6,74 

6.36 1 

**19 

2671.8 

-.09 

1.33 

0 

-.09 

1.33 ! 

6.32 

5.85 1 

*Ml(b) 

50.30a 

0 

0 

0 

0 

0 i 

10.66 

10.66 1 

*M2 

53.74a 

.66 

CO 

• 

.035 

-.07 

.84 1 

6.93 

6.61 ! 

*M3 

54.96a 

1.57 

1.20 

.10 

.11 

.96 1 

3.87 

3.83 1 

*M4 

57.95a 

.03 

1.33 

0 

.03 

1.33 1 

7.45 

6.93 1 


a - values in GHz 
b - used to determine eraissivity 

* - used in the GLAS retrieval systems for tenperature sounding 
** - 21 night cases only 
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,ccn 5 )uted brightness ten^Jeratures is due to sanpling differences betvreen radio- 
sonde and satellite observations in space and time and errors in the radio- 
sondes reports themselves. 

It IS apparent fran column 3 of Table 1.2 that significant bias errors 
exist in a number of the sounding channels. Biases in channels sensitive to 
the surface or upper stratospheric tsiperatures may result in part fran incor- 
rect specifications of the appropriate conditions. Nevertheless, it is clear 
that atmosEheric attenuation is underestimated in a nun±>er of channels, causing 
computed brightness temperatures to be systematically warm. 

As shown in previous studies [Jastrcw and Halem, 1973 ; ^feinreb, 1979] , 
systematic errors can be removed by empirical tuning of the transmittance func- 
tions. Vfe considers the following form 

T^(P) = T^(P)1+®1 (1.12) 

which is equivalent to multiplying the effective channel optical depths by 
1 +^ 1 * 

Tuning coefficients were found for channels, 4-7, 13-15, and M2-M4 vAiich 
manimized the bias in the brightness temperatures errors computed in a sample 
of 20 clear cases taken from the first two weeks of January. The tuning coef- 
ficients obtained for these channels are shown in column 5 of Table II. The 
mean and standard deviation of the errors of the computed brightness tempera- 
tures using the tuned transnattances for the vrtiole period are shown in columns 
6 and 7 of the table. No attempt was made to tune channels sensitive to the 
stratosphere or the surface. The empirical tuning coefficients found in 
Table 1.2 are similar to those found by Weinreb (1979). In all cases, the 
tuning coefficients are positive, indicating the need for additional attenua- 
tion to that included in the calculations. 
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The tuning coefficients determined from data of the first two we^s of 
the period leave small biases vrfien used on the whole sanple. The standard 
deviation of the errors remain basically unchanged except for channels 4 and M3 
viiich required large tuning coefficients. In these channels, significant 
uiprovements in the standard deviation occured as well. TSie tuning coefficient 
in channel 7 appears to be somewhat too large viien judged against the whole 
set. With the exception of the stratospheric sounding channel 2, the standard 

deviation of the tuned 15 vm and 4. Sum brightness tenperatures in channels used 
for retrievals of tenperature profiles are of the order of .7°. 'Qiis variahce 

in the agreeitent of computed and okserved brightness temperatures is comparable 
with the criterion used for the convergence of me physically derived tempera- 
ture profiles. Agreement in channel 7, used only for cloud height determina- 
tion, IS somewhat worse, possibly because of the effects of incorrect surface 
temperature. Vfe believe the agreement in the microwave channels, which is 
poorer than the infra-red channels by about 50%, is the result of poorer colo- 
cations with the lower spatial resolution microwave observations. The eighth 
column of Table 1.2 shows the standard deviations of the measurements for each 
channel. The ratio of the standard deviation of the observations to the stan- 
dard deviation of the errors, vAiich is significant for all but channel 1, may 
be thought of as a signal to noise ratio for each channel. The last column 
shows the standard deviations of the computed brightness temperatures using 
the tuned transmittances. The computed brightness temperatures show basically 
the same variances as the observations. 

Cloudy Cases 

Table 1.3 gives similar statistics for cases in the Jan-Feb 1979 period in 
which both fields of view were found to be partially cloud covered, that is, 

> 0. In this case, brightness temperatures computed from the radiosonde pro- 
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file vrere coonpared to brightness tenperature reconstructed from the satellite 
observations using equations (1.10) and (1.11). Tb avoid highly overcast situa- 
tion, only cases with n _< 2 were included in the statistics. Because n is 
determined by channel 13, radiances in this channel give "perfect agreement". 

No cloud correction is made to the microwave observations. As in Table 1.2, 
statistics for the window channels are given for the night cases only. 

All ccmputations were done using transmittance s tuned with the coeffi- 
cients determined from the subset of clear observations occuring in the first 
two wedcs of January. Ihese coefficients are shown in column 2. Ihe next two 
columns show the mean and standard deviation of the difference between ccm- 
puted and reconstructed brightness temperatures. The fifth column gives the 
root mean square of the difference between the reconstructed brightness tem- 
perature and the brightness temperature observed in field of view 1. For the 
channels sounding beneath the cloud layer, this difference is a measure of the 
cloud correction applied to the radiances. For channels sounding above the 
clouds, the difference between the observations in both fields of view is 
assumed to be noise. In these channels, noninally 1-3, the clear column 
radiance is taken as the average of the radiances in both fields of view, 
using an effective t of -.5. The difference between reconstructed and observed 
radiances is then a measure of instrumental noise in these channels. As in 
Table 1.2, the last two columns show the standard deviation of the reconstructed 
clear column and coitputed brightness tenperatures. 

The bias errors shown in Table 1.3 are all small and coiparable to those 
found in the clear cases. The tuning coefficients determined from clear cases 
are therefore applicable to cloudy cases as well. The standard deviation of 
the errors average .25® C larger than those in clear cases. The increase 
in error is small compared to the cloud conrections made, however, and indi- 
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TABIB 1.3 


COMPARISON OF RECONSTRUCTED AND COMPUTED BRIGHTNESS TEMPERATURES (°C) 
CLOUDY CASES JAN - FEB 1979 (146 CASES) 


CH 

e(a) 

MEAN 

(OOMP-RECON) 

STD DEV 

RMS 

REOON-OBS 

STD DEV 
REOON 

STD EGV 
COMP 

1 

0 

CO 

. 

2.60 

1.34 

3.99 

4.46 

*2 

0 

-.40 

1.31 

.38 

5.13 

5.06 

3 

0 

.91 

1.46 

.29 

5.02 

4.40 

*4 

.12 

-.02 

1.11 

1.11 

3.21 

3.27 

5 

.03 

.02 

1.38 

2.55 

4.83 

4.97 

6 

.015 

-.23 

1.64 

3.86 

5.89 

6.08 

7 

.04 

-.08 

1.85 

5.11 

5.83 

5.78 

**8 

0 

.42 

1.19 

5.76 

7.01 

6.68 

*13 (c) 

0 

0 

0 

3.90 

8.09 

8.09 

*14 

.04 

-.02 

.60 

3.17 

8.19 

8.50 

*15 

.035 

-.04 

1.05 

2.29 

6.83 

7.19 

16 

0 

.14 

2.51 

.65 

3.36 

4.03 

**18 

0 

.22 

1.18 

3.54 

9.44 

8.97 

**19 

0 

.23 

1.25 

3.33 

8.91 

8.39 

Ml(b) 

0 

0 

0 

0 

10.75 

10.75 

*M2 

.035 

-.07 

.85 

0 

7.60 

7.42 

*M3 

o 
1— 1 

• 

-.17 

.98 

0 

4.56 

4.30 

*M4 

0 

-.43 

1.43 

0 

8.50 

7.86 


a - e determined from 20 clear Jan cases 
b - used to determine e emissivity 
c - used to determine n 
* - used for temperature sounding 
** - 64 night cases only 
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cate that infra-red observations can be used in partially cloudy areas. The 
degradation in accuracy is scxtevrfiat worse in channels 5, 6, and 7 vdiich are 
not used in teirperature sounding. The agreanent in the window channels/ vrtiich 
IS in part determined by the accuracy of the climatological sea surface tenp- 
eratures, is not degraded under partially cloudy conditions. 

Tables 1.4 and 1,5 are equivalent to Tables 1.2 and 1.3 for cases taken from 
May 1979. The average absolute value of the bias in the untuned case is ,73* C 
and the average of the standard deviation is .96° C. In both sets, the tuning 
coefficients obtained from the clear cases- in the first two weeks of January 
were used in the calculations. The untuned biases shown in Table 1.4 are very 
similar to those in Table 1.2 for those channels not sensitive to the strato- 
sphere or surface tenperature. Consequently, the transmittances tuned with the 
set of winter tuning coefficients produce only small residual biases in the 
conputed brightness teitperatures. Tuning of channel 4 and M3 again decreased 
the standard deviation of the errors as well as removing the biases. As in the 
winter case, it appears that the tuning coefficient for channel 7 is somevrfiat 
too large. The standard deviations of the errors in the spring cases are simi- 
lar to those found in the winter. Channels 2, 3, and M4, sounding the mid-lower 
stratosphere, all appear to have a significantly colder bias in the spring 
than in the winter. This may be reflective of errors in the spring climatology 
used in extrapolating the radiosonde profiles above their highest reported 
level. The NMC stratospheric analysis was not available for use with the spring 
cases. The results in Table 1.5, showing statistics for cloudy cases from the 
spring period, are likewise very similar to those of Table 1.3. The average 
standard deviation of the errors is .3* higher in the May cloudy cases than in 
the May clear cases, but as in the winter, the largest increases in error occur 
in channels not used in the GLAS tonperature retrievals. As in Table 1.4, the 
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TAHLE 1.4 


COMPARISON OF OBSERVED AND OC»!E>UTED BRIGHTNESS TEMPERATURE (®C) 
CLEAR CASES MAX 5 - MAZ 20 1979 (33 CASES) 

UNTUNED TUNED 


CH 

1 

MEAN 

COMP-OBS 

STD c®:v 

e(a) 

MEAN 

STD DEV 

STD I®V 
CBS 

STD DEV 
GCMP 

1 

1 1 

1.55 

2.00 

0 

1.55 

2.00 

2.57 

2.72 

1 

*2 j 

-1.55 

1.14 

0 

-1.55 

1.14 

2.97 

3.10 

1 

3 1 

t 

.64 

1.00 

0 

.64 

1.00 

2.46 

2.36 

1 

*4 1 

1.62 

.72 

.12 

-.27 

.67 

2.20 

2.32 

1 

5 j 
1 

.54 

.67 

.03 

-.10 

.67 

3.58 

3.61 

1 

6 1 

.18 

.81 

.015 

-.25 

.81 

4.24 

4.28 

1 

7 1 
1 

.32 

1.02 

.04 

-.37 

1.02 

4.14 

3.96 

**8 i 
1 

-.54 

1.18 

0 

-.54 

1.18 

3.83 

3.54 

1 

13 1 

-.14 

.92 

0 

-.14 

.92 

5.67 

5.33 

1 

14 1 

.72 

.66 

.04 

-.09 

.67 

5.79 

5.70 

1 

15 1 

.80 

.62 

.035 

-.02 

.64 

4.98 

4.95 

! 

16 1 
1 

.87 

1.84 

0 

.87 

1.84 

1.97 

2.67 

**18 i 

I 

-.39 

.56 

0 

-.09 

.56 

4.87 

4.81 

**19 1 

-.38 

.61 

0 

-.02 

.61 

4.50 

4.45 

1 

Ml 1 

0 

o(t>) 

0 

0 

0 

8.91 

8.91 

1 

*M2 1 

.30 

.72 

.035 

-.44 

.71 

5.89 

5.69 

1 

*M3 1 

1 

1.32 

.91 

.10 

-.22 

.79 

3.65 

3.66 

I 

*M4 1 

-.61 

.97 

0 

-.61 

.97 

4.53 

4.39 


a - ^ determined frcsii 20 clear Jan. cases 
b - used to determine omissivity 
* -• used for tenperature sounding 
** - 19 night cases only 
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TABLE 1.5 


COMPARISCW OF RECmSTRUCTED AND COMPUTED BRIGHTNESS TEMPERATURES (°C) 
CLOUDy CASES MA^ 5 - MKiT 20 1979 (47 CASES) 


CH 

e (a) 

MEAN 

(CCMEMJBS) 

STD DEV 

RMS 

RECON-OBS 

STD DEV 
REOON 

STD DEV 
COMP 


0 

1.66 

2.08 

1.42 

3.68 

3.42 

*2 

0 

-.84 

1.28 

.38 

4.01 

4.16 

3 

0 

.30 

1.38 

.30 

4.35 

3.69 

*4 

.12 

-.09 

1.04 

.73 

1.96 

1.77 

5 

.03 

.09 

1.50 

1.67 

3.81 

3.81 

6 

.015 

-.31 

1.73 

2.63 

5.31 

5.43 

7 

.04 

-.32 

1.84 

3.79 

5.43 

5.63 

**8 

0 

.31 

.83 

5.18 

6.12 

6.24 

*13 (c) 

0 

0 

0 

2.01 

7.60 

7.60 

*14 

.04 

.01 

.51 

1.37 

7.42 

7.53 

*15 

.035 

.34 

1.05 

.43 

5.83 

5,76 

16 

0 

-.57 

1.81 

.08 

2.22 

1.80 

**18 

0 

-.35 

.72 

3.42 

7.64 

7.57 

**19 

0 

-.72 

1.31 

3.29 

7.42 

7.53 

Ml(t>) 

0 

0 

0 

0 

11.53 

11.53 

*M2 

.035 

.06 

.70 

0 

6.10 

5.89 

*M3 

.10 

-.36 

.84 

0 

3.80 

3.51 

*M4 

0 

-.83 

1.23 

0 

6.82 

6.39 


a - e determined from 20 clear Jan. cases 
b - used to determine emissivity 
c - used to determine t 
* - used for tenperature sounding 
** - 18 night cases only 
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biases for channels 2 , 3/ and M4 are again significantly lower than in the 
winter period. 

1.7 . SiMnary 

Tables 1.2-1. 5 show that under clear ccaiditions, given the tenperature-huna- 
dity profile, radiances can be calciiLated which agree with observations to 
about .7°C for the HIRS2 channels not sensitive to the upper stratosphere or 
the surface, and about 1°C for the MSU channels. Part of these differences 
are due to noise and sampling differences in space and tune between the satel- 
lite and radiosonde. The accuracy of this calculation is sufficient to allow 
for determination of atmospheric temperature profiles via an iterative physical 
technique m which teirperature profiles are found, which when substituted in 
the radiative transfer, match the observations to a desired accuracy. Under 
partially cloudy conditions, cloud corrections can be made to the infra-red 
observations to still allow for accurate temperature soundings. 

Significant bias errors exist in the calculations for soite channels, 
however. These errors can be ranoved to a large extent by an eitpirical tuning 
of the transmittance functions with coefficients vrfiich appear to be transfer- 
able fron one season to another. In all cases, additional atmospheric attenua- 
tion appears necessary in the calculated transmittance functions. 

While empirical tuning coefficients are sufficient for practical use, it 
is much more desirable to directly eliminate that portion of the bias that 
results from theoretical uncertainties. Coffee and Goldman (1981) have indi- 
cated that atmospheric absorption due to the V 2 band of HCN, centered at 
712 cm"*l, IS significant in that spectral region. Absorption by HCN was not 
included in the line by line calculations done in this work. In addition, a 
more accurate value of .31 ppmv for the tropospheric mixing ratio of N 2 O 
[Chedin, 1982] , rather than the value of .28 ppmv, reported in McClatchey ^ 
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al « [1973] should be used in the calculations. Hfe are currently including HCN 
absorpticai in our calculation as vrell as modifying the N 2 O mixing ratio to 
see the extent to vrtiich these factors may explain the untuned biases in HIRS 
channels S-'S, and 14-15. 

The rapid transmittance algorithm used in the calculations in described in 
the next chapter. It will be shown to contribute only small errors to the radiance 
calculations as cotpared to calculations based on line-by-line ccmputed trans- 
mittance functions. 
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Chapter 2 


A Rapid Algorithm for Modeling Atmospheric Transmittance s 
2.1 Tlie physical basis for the form of the rapid algorithm 

Hie rapid transmittance algorithm is of the form of equation 1,7, Hie 
averaged discrete line transmittance through the atmosjAiere from pressure 
Pjj to the top of the atmosphere, at a zenith 0 , as seen by channel i, is 
modelled as 


if(^ 3 ^^3“ 1^® ) '^i0t^3 f^3“l^® ) ^iw(^3 ^^3— 1»^® 

3=1 


( 2 , 1 ) 


where and represent models for effective layer transraittances 

from pressure to P 3 -I (P^>P^..l) at zenith angle 0, Hie term t” 3 ^p represents 
absorption by gases assuned to have a fixed mixing ratio, while ^ ^iW 
represent absorption due to ozone and water vapor respectively, ^iL(Pjlf®) 
from equation ( 2 , 1 ) is used to model ■riij(P, 0 ) from equation 1 , 2 , 

Line-by-line calculations done at zenith angles of 0®, 50®, and 70°, for 
eight climatological profiles shown in Table 2, 1-2. 3 are used to generate the 
coefficients for the effective transmittance models at the appropriate angle. 
Effective layer transmittances at other zenith angles are obtained by linear 
interpolation of the logarithm of the effective layer transmittance in sec 0 . 
Til(Pjj_, 0 ) coTputed frcm equation 2.1 is then multiplied by additional factors 
to account for H 2 O continuum absorpton, N 2 continuum absorption, and scattering 
as discussed in the text. 

Because channel 1 is not monochromatic, the effective layer transmit- 
tances do not obey the properties of equation 1.4. Instead, given line by 
line transmittance calculations for T 2 p(P, 9 ) , tiEDt^'®)' ^ 

correspondirig respectively to absorption using only gases of fixed distribution, 
using fixed gases and ozone, and using all species, we define effective mean 
layer transmittances 
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TABLE 2.1 


CLIMATOLOGICAL TEMPERATURE PROFILES (®C) 


c 

•V 

P(MB) 

i 1 

1 JAN 0®N 1 

! 

JAN 20°N 1 

1 

JAN 40 °N 1 

JAN 50 ®N i 

JAN 70 “N 1 

JUL 20°N 1 

JUL 40°N 1 

1 

JUL 60°N 1 

STD. PRO 

«> 

1 

1 r 

I 270.0 I 

t 

270.0 1 

1 

1 

265.0 1 

1 

250.0 1 

235.0 1 

270.0 1 

285.0 1 

r 

290.0 1 

1 

265.0 


2 

I 263.4 I 

I 

263.7 1 

1 

259.0 1 

1 

245.8 1 

231.2 1 

264.1 1 

277.5 1 

282.4 1 

1 

259.7 


5 

I 241.4 1 

1 

242.7 1 

1 

239.0 1 

1 

231.8 ! 

218.2 1 

244.6 1 

252.5 1 

257.3 1 

1 

242.2 

A. 

10 

1 227.0 j 

1 

229.0 1 

I 

225.9 1 

1 

222.6 ! 

210.6 1 

231.9 1 

236.1 I 

240.9 1 

230.8 

v«t- 

20 

1 222.0 1 
1 1 

1 

222.5 1 

1 

219.7 1 

I 

217.7 1 

206.5 1 

226.5 t 

229.5 1 

1 

232.5 1 

1 

224.5 


30 

1 217.5 1 

1 

219.5 1 

1 

217.5 i 
1 

217.2 1 

206.0 ! 

222.5 1 

222.5 1 

1 

229.5 1 

221.0 

to 

00^^ 

50 

I 207.0 1 

1 

211.1 1 

214.9 ! 

j 

217.0 1 

206.8 1 

213.7 1 

219.1 i 

1 

226.8 1 
1 

213.5 


70 

i 196.6 1 

201.3 ! 

I 

214.1 1 

I 

217.4 j 

208.5 1 

206.6 1 

214.9 1 

1 

226.1 ! 
1 

208.1 


100 

1 193.7 1 

j 

198.1 1 

1 

214.0 1 

I 

218.8 1 

211.6 ! 

200.3 1 

212.3 1 

1 

225.6 1 

1 

206.1 


150 

i 204.0 1 

1 1 

208.1 1 
1 

217.5 ! 

I 

220.2 1 

213.4 1 

207.1 1 

215.6 1 

226.4 1 

1 

211.0 


200 

! 217.8 ! 

1 

[ 

218.1 1 
1 

217.9 1 

1 

218.2 1 

214.7 1 

218.9 1 

221.5 1 

1 

224.4 1 

1 

216.9 


250 

t 234.3 ! 

1 1 

227.4 1 

1 

1 

221.7 i 

j 

220.7 1 

215.0 1 

229.7 ! 

229.2 1 

1 

227.4 1 

1 

224.6 


300 

1 240.0 1 

1 1 

236.4 1 

1 

225.8 1 

1 

221.6 ! 

215.2 j 

239.9 i 

237.8 1 

230.8 1 

1 

232.1 


400 

1 255.3 1 

1 

251.3 1 

1 

238.1 1 

1 

232.7 ! 

225.9 1 

245.9 1 

252.0 1 

244.0 ! 

1 

247.3 


500 

1 266.7 1 

1 

262.7 1 

1 

1 

248.1 1 

1 

241.9 ! 

235.0 1 

266.2 1 

262.9 1 

254.8 1 

1 

259.0 


700 

1 282.2 1 
j 1 

279.1 1 

1 

264.6 1 

\ 

247.6 1 

249.2 ! 

282.9 1 

279.5 1 

270.7 1 

1 

275.9 


850 

1 289.5 1 

j 1 

285.5 ! 

1 

1 

270.8 1 

j 

262.8 1 

255.4 1 

291.6 1 

288.3 I 

279.3 1 

284.2 


1000 

1 295.5 i 

i 1 

292.5 1 

1 

277.7 1 

1 

270.1 1 

265.8 1 

297.8 i 

292.0 1 

I 

284.5 1 

1 

290.0 



TABLE 2.2 


CLIMATOLOGICAL OZC»3E PROFILES (PPMV) 


P(MB) 

i 1 

1 JAN 0°N 1 
1 P 1 

JAN 20®N 
P 

1 

1 

1 

1 

JAN 40 °N 1 
P 1 

1 

JAN 50°N 1 
p ! 

i 

JAN 70®N 1 
P 1 

1 

JUL 20°N i 
P 1 

JUL 40®N 1 
P 1 

JUL 60°N 
P 

1 

1 r 

1 3.818 1 

1 

3.812 

T 

1 

2.661 i 
1 

r 

2.023 1 

t 

1.342 1 

1 

r 

3.805 ! 

1 

2.352 1 

1.547 

2 

1 3.818 i 

3.812 

1 

1 

2.661 1 

1 

2.023 1 

1 

1.342 1 

1 

3.805 I 
1 

2.352 1 

1.547 

5 

\ \ 

1 4.583 1 

1 1 

4.558 

1 

3.395 ! 

I 

2.657 1 

1 

1.744 1 

4.348 1 

1 

3.067 1 

2.020 

10 

1 1 

1 6.876 1 

1 1 

6.790 

1 

5.599 1 

1 

4.559 1 

1 

2.952 ! 

5.978 1 

1 

5.211 1 

3.439 

20 

1 1 
1 7.808 1 

1 

7.778 

1 

1 

6.903 1 

1 

5.517 1 

1 

3.889 1 

1 

6.542 1 

1 

6.180 i 

3.979 

30 

1 7.119 t 

1 1 

7.209 

1 

6.446 1 

1 

1 

5.496 ! 

1 

4.084 1 

1 

1 

5.978 ! 

1 

5.322 j 

3.866 

50 

1 2.989 1 

1 1 

3.774 

1 

4.438 1 

1 

4.711 i 
1 

4.318 1 

1 

3.080 ! 

1 

3.291 1 

3.140 

70 

1 1 
1 .634 1 

1 

1.449 

1 

2.446 1 

j 

3.351 1 

3.623 1 

1 

1.419 1 

1 

2.023 i 

2.597 

100 

t .242 1 

1 1 

.513 

1 

1.147 1 

j 

2.144 1 

1 

2.838 1 

1 

.604 1 

1 

1.087 1 

1.932 

150 

I .085 1 

1 1 

.154 

1 

.664 ! 

1.208 1 
1 

1.902 1 

1 

.205 1 

1 

.483 1 

1.177 

200 

1 1 
j .042 1 

1 

.063 

1 

1 

.332 1 

.694 1 

1 

1.178 1 

.094 1 

1 

.272 i 

.694 

250 

1 .034 1 

1 1 

.053 

1 

.168 I 

1 

1 

.400 ! 

1 

.641 ! 

1 

.071 1 

.193 1 

.398 

300 

1 .031 1 

i t 

.046 

! 

.145 1 

1 

.254 1 

1 

.374 1 

1 

.064 i 

.145 1 

.253 

400 

1 1 
1 .030 1 

1 

.038 

! 

.101 1 
1 

.145 ! 

1 

.189 1 

1 

.057 1 

1 

.104 1 

.151 

500 

1 .030 1 

1 1 

.029 

1 

1 

.062 i 
1 

.072 1 

1 

1 

.083 ! 

1 

.051 1 

.071 1 

.091 

700 

1 .021 1 

.020 

I 

.033 1 

1 

.032 1 

1 

.026 ! 
1 

.039 1 

1 

.050 1 

.047 

850 

1 1 
1 .018 1 

.018 

1 

1 

.025 1 

1 

1 

.025 1 

.024 1 

.031 1 
1 

.048 1 

.037 

1000 

1 .018 1 
1 I 

.018 

1 

JL 

.024 1 

1 

.024 1 

1 

.024 i 
! 

1 

.026 1 
1 

.044 1 

.036 

Total O3 .244 

.282 


.355 

.430 

.480 

.263 

.307 

.362 


(cm atm) 



TABLE 2.3 


CLIMATOLOGICAL WATER VAPOR PROFILES (gm/cm2atm) 


P{MB) 

I JAN 0°N 

JAN 20®N 

i JAN 40®N 

JAN 50 °N 

JAN 70°N 

1 JUL 20®N 

1 JUL 40°N 

i JUL 60®N 

1 


1 P 

P 

1 P 

P 

P 

1 P 

1 P 

i P 

J 

1 

1 

1 .002 

.002 

1 .002 

.002 

.002 

1 .002 

1 .002 

1 .002 

1 

1 

2 

1 .002 

.002 

1 .002 

.002 

.002 

1 .002 

I .002 

i .002 

1 

1 

1 

5 

! .002 
1 

.002 

1 .002 

.002 

.002 

1 .002 

1 .002 

1 .002 

1 

1 

10 

1 .002 
1 

.002 

1 .002 

.002 

.002 

1 .002 

1 .002 

1 .002 

1 

20 

1 .002 

.002 

1 .002 

.002 

.002 

1 .002 

1 .002 

1 .002 

1 

1 

30 

i .002 

.002 

1 .002 

.002 

.002 

1 .002 

1 .002 

1 .002 

1 

1 

1 

50 

1 .002 

.002 

1 .002 

.002 

.003 

! .002 

1 .002 

1 .002 

1 

1 

1 

70 

1 .002 
1 

.002 

1 .002 

.003 

.006 

! .002 

1 .002 

1 .002 

1 

1 

i 

100 

1 

1 .002 
1 

.002 

1 .002 

.007 

.010 

! .002 

1 .002 

1 .002 

1 

I 

1 

150 

1 

1 .006 
1 

.005 

1 .005 

.016 

.017 

I .006 

1 .006 

1 .002 

\ 

1 

1 

200 

1 .027 

1 

.018 

1 .013 

.032 

.029 

! .028 

1 .023 

1 .002 

1 

1 

1 

250 

1 .11 

.06 

I .05 

.07 

.04 

1 .12 

1 .09 

1 .01 

1 

! 

1 

300 

1 .21 
j 

.11 

1 .06 

.09 

.06 

1 .23 

! .16 

1 .03 

1 

1 

1 

400 

1 .59 

1 

.25 

1 .13 

.15 

.08 

1 .65 

1 .43 

1 .18 

1 

1 

1 

500 

1 1.48 

1 

.66 

1 .36 

.26 

.15 

1 1.68 

1 1.12 

1 .71 

1 

1 

1 

700 

1 4.01 

1 

2.53 

! 1.28 

.84 

.43 

1 4.92 

1 3.60 

1 2.55 

1 

1 

1 

850 

1 8.27 

1 

5.95 

1 2.63 

1.70 

.76 

1 9.81 

1 7.09 

1 4.85 

1 

1 

1 

1000 

i 12.38 
! 

9.25 

! 3.95 

2.31 

.82 

1 14.42 

1 10.72 

1 6.84 

1 

1 

1 

Total HoO 3.64 
gm/cm2 

2.53 

1.14 

0.73 

0.32 

4.31 

3.14 

2.08 






( 2 . 2 ) 


^ 10(P: iFO(P:3,Pj-l,9 ) (2.3) 

•^lf(P3 ^3-1,9) 
and 

iw(P: ,^3-1,0 ) = iTOW(P3 ,^3-1,9 ) (2.4) 

iED(P3 ,^3“1,9 ) 

Equation 2.2 defines an effective mean layer transmittance based on line-l^- 
line calculations using any combination of constituents [Halem and Susskind, 

1977]. Note that the effective layer transmittance s for ozone (equation A3) and 
water vapor (equation 2.4) do not involve in any way calculations based on absorp- 
tion of water vapor or ozone alone [Susskind, 1975] . In fact, significant 
differences may exist between effective layer transmittances defined from 
equations 2.3, 2.4 and those defined in equation 2.2 based on the single species 
transmittances. Table A2 of Halem and Susskind [1977] shows that the bright- 
ness tenperatures conputed for VTPR channel 7 using line-ty-line transmittances, 
'^FO'^VIr differ from those computed using line-by-line '^fcW/ t>y .4‘>C for a 
tropical temperature humidity profile. Ihe spectral response of VTPR channel 
7 IS very similar to that of HIRS channel 7. Errors of similar magnitude are 
then expected for HIRS2 channel 7. 

2.2 The model for water vapor and ozone absorption 

"Uie basic assumption of the models of water vapor and ozone transmittance 
IS that the effective mean layer transmittances in eqs. (2.3, 2.4) can be treated 
as having the transmittance properties of a gas in a homogeneous layer having 
the mean teirperature T, and pressure P, of the atmospheric layer, and vertical 
column density u of the absorbing gas in the layer. This assumption is reason- 
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ably valid because use of equation 2,2 removes most of the dependence of the 
mean layer transmittance cn the properties of the atmosphere above the layer, 
and absorption due to water vapor and ozone has a secc«id order effect on the 
radiances in the tenperature sounding channels. We then expect the log of the 
mean layer transmittance to be pr<^»rtional to u for weakly absorbing lines, 
and uV2 for strong lines. For a conposite of lines, an effective exponent 
of intermediate value is obtained. Ihe absorption coefficient depends on the 
pressure P and the temperature T. 

The following form was therefore used to model the effective water and ozone 
transmittances for all channels and all layers: 

^ic^^3^^3-l,®^ " tl - Bj^^^,(T^-273)] Uj,(3,3-l)^i,c (2.5) 

where c stands for constituent, either ozone or water vapor, Uc(3^3“l) 
integrated colimin density of the species in the layer between 3 and 3-I, 
is a channel and species dependent constant betvreen .5 and 1, Ai^-j^c an 
effective channel, species, pressure, and angle dependent absorption coeffi- 
cient, and ^ channel and^^cies dependent constant (pei>cent change 

per degree). For simplicity, the temperature dependence, ^<3 exponent, 

^x,c' taken to be independent of pressure and angle. The coefficients A, 

B, and N are determined from the effective mean layer transmittances computed 
from the line-by-line calculations from the eight climatological profiles 
shown in Tables 2.1— 2.3. Table 2.4 shows select values of the coefficients for 
the water vapor and ozone models. A(P,6) is raised to the power 1/N so as to 
be in units of (mol/cm''^)*"^. The absorption coefficients are shown at 1000 mb 
and 100 mb respectively, \diere absorption due to water vapor and ozone is most 
significant. Because of the complex nature of the effective transmittance, the 
pressure behavior of A is not simple. The range of coefficients for a parti- 
cular channel in the pressure region v^ere significant absorption is present 
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TABIiE 2.4 


SELECT WATER VAPOR AND OZC»ffi TRANSMITTANCE COEFFICIENTS 

WATER VAPOR OZCNE 


CHANNEL 

A(1000.0)1/N 

{M0L-1CM2) 

B 

N 

A(100,0)1AI 

(MOL“1cm2) 

B 

1 

0 

0 

1.0 

0 

0 

2 

0 

0 

1.0 

6.2 E-23 

0 

3 

0 

0 

1.0 

2.6 E-21 

0 

4 

1.6 E-24 

0 

.8 

3.6 E-22 

0 

5 

4.2 E-25 

0 

.85 

1.9 E-21 

0 

6 

8.4 E-24 

0 

.65 

5.1 E-21 

0 

7 

4.4 E-26 

.016 

.5 

3.1 E-21 

0 

8 

4.5 B-25 

0 

.9 

0 

0 

9 

6.9 E-25 

0 

1. 

7.8 E-20 

0 

10 

1.1 E-24 

0 

.8 

0 

0 

11 

5.1 E-23 

0 

.7 

0 

0 

12 

4.2 E-22 

0 

.6 

0 

0 

13 

1.4 E-24 

0 

1.0 

0 

0 

14 

2.8 E-25 

.019 

.8 

0 

0 

15 

6.8 E-25 

.017 

.95 

0 

0 

16 

1.4 E-24 

0 

1.0 

0 

0 

18 

6.7 E-26 

.016 

1.0 

0 

0 

19 

3.9 E-25 

0 

.8 

0 

0 

Ml 

1.1 E-24 

0 

1.0 

0 

0 

M2 

1.2 E-24 

0 

1.0 

0 

0 

M3 

1.3 E-24 

0 

1.0 

0 

0 

M4 

1.4 E-24 

0 

1.0 

0 

0 


N 

1 

.7 

.8 

1 

1 

1 

.95 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
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IS generally less than 20%, Por most channels, B is found to be nearly zero 
and N IS nearly 1, 

McMillin et al ., (1979), and Weinreb and Neuendorfer (1973) have given 
more elaborate models to account for the effects of water vapor on atmospheric 
transmittances. In both cases, the coefficients for the models are based on 
calculations of averaged transmittances due to water vapor alone, Uiese modeled 
transmittances are multiplied by the dry transmittances to give total trans- 
mittances as in e(^ation 2,1, The appropriate quantities bo use are the effec- 
tive mean layer transmittances as defined in equation 2,4, v^ich differs signi- 
ficantly from the mean layer transmittances ccHtputed from the absorbing gas 
alone, because absorption by a gas such as H 2 O is not totally uncorrelated 
with absorption by the fixed gases, especially in a relatively narrow channel, 
where only one or two H 2 O lines may be of significance, 

2.3 The model for absorption due to gases of fixed distribution 

Most of the absorption for the teirperature sounding channels is due to the 
gases of fixed distribution, primarily OO 2 and N 2 O. Hie transmittance at a 
given angle depends only on the tenperature profile. Ohe effective mean layer 
transmittance for each reference angle is modelled according to 

V - 11,°) + Cy(e)(T^(e)^^^0(e) ^^. 6 ) 

vAiere is the mean tenperature in the layer between Pj to P^-i, for the 
tenperature profile under consideration, is the mean tenperature in layer 
:) in a standard teirperature profile shown in Table Al, and T and T ° are 
effective mean teirperatures for the entire profile froti P-^ to the top of the 
atmosphere for the tenperature profile under consideration and the standard 
temperature profile respectively. The effective mean temperature above pres- 
sure p for channel i is defined as the average tarperature above pressure P 
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weighted by the weighting function for channel i. The effective tenperature 
IS then channel and angle dependent and is defined as 

P 

Tn(®) = 1 / \(P) dTT{9) dP (2.7) 

l-f.^(P^ w) o dP 

where (Pf®) is the transmittance of channel i for the standard teirperature 
profile. 

The coefficients ( 6 ) / ( 9 ) , and ( 6 ) are determined so as to give 

the best fit in the least squares sense to the values of iip obtained fron 
line-by-line calculations. As expected , the coefficient Ai^(0) was found to 
be very dose to 'Cjp°(P^^P^_ 2 ^^ 0 ) , the effective layer transmittances for 
the standard profile shown in Table 2.1. 

The algorithm is essentially equivalent to that of McMillin and Fleming 
(1978), except that McMillin and Fleming ejqsand the effective layer transmit- 
tances in a four term expansion about a standard profile effective layer trans- 
mittance. Two of their terms deal with the tenperature in the layer and the 
remaining two contain channel independent vieighted temperatures above the 
layer. 

2.4 Ability of the model to fit the climatological profiles 

The inportance of the transmittance error is the error made in computing 
brightness temperature, or equivalent black bo<^ temperature, which arises fron 
using incorrect transirattances in equation 1 of the text. Table 2.5 shews the 
ability of the rapid algorithm to reproduce the brightness temperatures for the 
profiles used to generate the coefficients. The second and third columns of 
Table 2.5 summarize the maximum and RMS errors in brightness temperatures for 
the eight profiles computed using the transmittance functions from equation 2.6 and 
the coefficients determined from the fit as compared to those computed using 
F, the line-by-line transmittance. In ccmputing these brightness tempera- 
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TABLE 2.5 


TRANSMITTANCE MODEL BRIGHTNESS TEMPERATURE ERRORS ( °C) FOR HIRS CHANNELS AT NADIR VIEWING 



MOI®L ERRORS 

MOEEL 

ERRORS 

MODEL ERRORS 


MAXIMUM EFFECT 

CHANNEL 

FIXED GAS 

FIXED+OZCNE 

FIXEDK)ZC»^E+WATER 

azom 

1 WATER 

0ZQNE4WATER 


MAX 

RMS 

MAX 

RMS 

MAX 

RMS 



4<X)NTINUUM 

1 

.038 

.021 

.038 

.021 

.038 

.021 

0 

0 

0 

2 

.085 

.055 

.085 

.055 

.085 

.055 

0 

0 

0 

3 

.015 

.009 

.013 

.006 

.013 

.006 

.27 

0 

.27 

4 

.071 

.041 

.079 

.048 

.077 

.050 

.81 

.47 

1.13 

5 

.144 

.074 

.129 

.079 

.145 

.082 

2.04 

1.08 

2.90 

6 

.092 

.049 

.076 

.047 

.105 

.063 

1.58 

1.92 

4.57 

7 

.153 

.086 

.119 

.077 

.115 

.079 

1.10 

3.33 

5.80 

8 

.002 

.001 

.002 

.001 

.103 

.063 

0 

1.33 

4.89 

9 

.008 

.003 

1.535 

.922 

1.464 

.891 

26.81 

.70 

27.48 

10 

.003 

.002 

.003 

.002 

.264 

.122 

0 

6.26 

7.71 

11 

.005 

.003 

.005 

.003 

1.108 

.605 

0 

33.23 

33.23 

12 

.000 

.000 

.000 

.000 

1.793 

.982 

0 

57.57 

57.57 

13 

.014 

.007 

.014 

.007 

.073 

.038 

0 

1.14 

1.18 

14 

.027 

.012 

.027 

.012 

.028 

.014 

0 

.57 

.58 

15 

.033 

.014 

,033 

.014 

.033 

.013 

0 

.30 

.31 

16 

.085 

.039 

.085 

.039 

.085 

.041 

0 

.06 

.06 

18 

.003 

.001 

.003 

.001 

.003 

.001 

0 

.08 

.10 

19 

.001 

.000 

.001 

.000 

.050 

.039 

0 

1.71 

1.73 


ORIGINAL PAGH IS 
OF POOR QUALITY 


tures. It was assumed that the surface was a black body whose tejtperature was 
equal to the 1000 mb air tenperature. Ihe maximum difference in brightness 
temperature arising from use of the modelled transmittances rather than the 
exact transmittances is equal to .15®C and is on the order of the instrumental 
noise. The next four columns indicate maximum and RMS errors for the bright- 
ness teiTperatures coirputed for the eight profiles using ipo c^pared to 
Tp tq and ■'^FCW canpared to tp Xq ’’■W* exact tenperature-humidity-ozone 
profiles were used in the calculations. 

The last three columns show the extent to which absorption due to O3 and 
water vapor effects the brightness tei^ratures of each channel. The column 
marked maximum O3 effect gives the maximum difference of brightness tenperature, 
conputed/ on the one hand, using line by line transmittances calculated with 
only the fixed gases, Xp, and, on the other hand, using line-l^-line trans- 
itu, ttances calculated also includir^ O3 absorption, x^^, column marked 

maximum H2O effect gives the brightness temperature difference ccnputed using 
xpo on the one hand, and tpcw, the transmittances including all species, 
on the other. The maximum total effect, including absorption due to ozone and 
water vapor lines and water continuum, is shown in the last column. 

Except for channels 8-12, which are sensitive primarily to absorption of 
H2O or O3 themselves, the transmittance model fits the eight climatological pro- 
files extremely well, the largest source of error coming from the fixed gas mo- 
del. In channels, 9, 11, and 12, where the effects of absorption due to O3 and 
H2O are the order of 20-55°C, the rapid algorithm errors are larger than 1®C and 
a more elaborate algorithm appears to be necessary to reduce the errors. Inclu- 
sion of a term representative of the effective species column density above the 
layer, analogous to the third term in equation 2.6, as well as allowing for a 
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pressure dependence of N, should achieve more accurate results. 

2,5 Test of the model on independent profiles 

Line-by-line calculations vjere run for 2 of the 72 clear profiles from the 
January - February period, to assess the contribution of the rapid transmittance 
algorithm as a source of error in the radiance calculations for real profiles. 
The first was a mid-latitude profile chosen because the tenperature profile 
differed significantly frcm the climatological set used to generate the rapid 
algorithm coefficients. Tlie second profile was a very humid trcpical profile 
observed at a large zenith angle. Hi is was chosen to test the water vapor cor- 
rection model in an extr^ne case, Hie two profiles are shown in Table 2,6, The 
radiosonde temperature profiles were reported to 10 mb and 20 mb respectively. 
Hie humidity profiles were reported to 400 and 300 mb. Values extrapolated from 
the radiosonde reports are marked by *. Climatological values, marked by **, 
were used for the sea surface tenperature, and the O 3 profile, the 1 mb tem- 
perature tenperature, and humidities above 100 mb. 

Table 2.7 shows for each profile, Togg, the observed brightness taiperar 
tures for each channel, the difference between Tqbs and Tj^, the brightness 
tenperature s coirputed using equation 1.1 and the untuned rapid algorithm model 
generated transmittance coefficients, and the differences between and Tg, 
the brightness tenperatures conputed using line-by-line generated transmit- 
tances. Hie calculations are all run at the appropriate zenith angles, 22.8® 
for the mid-latitude profile and 46,9® for the tropical profile. Hius, all 
aspects of the rapid algorithm, including zenith angle interpolation, are 
tested simultaneously by this conparison. 

Hie brightness tenperature errors introduced by the rapid algorithm, shown 
in the columns Tj^Tg, are of the order of ,1®C for all channels but those prima- 
rily sensitive to water vapor or ozone themselves. The errors are similar to 
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TABLE 2.6 


ATMOSPHERIC SOUNDINGS USED TO TEST THE RAPID TRANSMITTANCE ADSORITHM 



MID LATITUIB PROFILE 


TROPICAL PROFILE 



FEB 9 OZ 47°N, 17°W 


JAN 5 OZ 14®S, 171®W 



ZENITH ANGLE = 22.8® 


ZENITH ANGLE = 46.9® 



SURFACE PRESSURE = 992 MB 

= 1.66 giTv'cm^ 

SURFACE PRESSURE = 1009 MB 



SST = **284.7®K PREC. WAT. 

SST = **302. 1®K PREC. WAT. 

= 6.56 gn\/cm 


O 3 PROFILES = **JAN 50®N 


OZONE PROFILE = **JAN 0°N 

P(MB) 

T(®K) 

g(gm/cm 2 atin) 

T®(K) 

q(gm/cm 2 atnt) 

1 

**259.1 

**.002 

**274.2 

**.002 

2 

*245.6 

**.002 

*258.2 

**.002 

5 

*227.6 

**.002 

*240.6 

**.002 

10 

208.7 

**.002 

*232.9 

**.002 

20 

206.1 

**.002 

225.7 

**.002 

30 

209.7 

**.002 

221.7 

**.002 

50 

215.1 

**.002 

205.9 

**.002 

70 

218.5 

**.002 

200.5 

**.002 

100 

222.1 

**.002 

191.7 

**.002 

150 

222.3 

*.005 

207.7 

*.02 

200 

224.7 

*.02 

221.5 

*.09 

250 

222.9 

*.04 

234.3 

*.28 

300 

218.1 

*.08 

244.5 

.76 

400 

235.5 

.24 

259.5 

2.0 

500 

248.1 

.85 

269.7 

3.9 

700 

267.3 

2.5 

283.2 

8.2 

850 

277.0 

4.7 

291.4 

16.3 

SURFACE 

285.4 

* 6.6 

300.2 

21.0 


**CLIMATOLOGy 

*EXTRAPOLATED 



TABLE 2.7 


ASSESSMENT OF ERRORS CAUSED BY RAPID ALGORITHM (°C) 

MID lATITUEG TROPICAL PROFILE 

PROFILE 


CHANNEL 

Tobs* 

“ Tqbs 


Tqbs 

'Ik - Tqbs 


1 

225.74 

-.16 

-.03 

234.28 

2.48 

-.10 

2 

217.48 

.43 

.06 

218.15 

1.21 

-.13 

3 

219.22 

.59 

.08 

216.69 

2.95 

.10 

4 

226.29 

1.28 

-.04 

227.56 

3.32 

.03 

5 

238.63 

.71 

-.08 

243.85 

1.25 

.02 

6 

250.23 

.13 

-.12 

255.02 

.90 

-.02 

7 

262.25 

.51 

-.02 

265.89 

1.20 

.07 

8 

282.75 

-.66 

-.10 

289.85 

-.56 

-.68 

9 

251.79 

2.44 

-.85 

272.53 

-2.96 

-1.43 

10 

278.81 

-.75 

-.03 

284.67 

-1.20 

-.81 

11 

255.86 

-5.64 

-.20 

253.64 

-1.53 

-.99 

12 

242.30 

-10.30 

.02 

241.31 

-5.96 

.50 

13 

268.49 

-.30 

.02 

278.47 

-.11 

-.12 

14 

256.56 

1.03 

.05 

265.67 

1.38 

.10 

15 

243.84 

1.17 

.08 

251.82 

1.41 

.06 

16 

225.82 

1.25 

.10 

231.44 

1.58 

.16 

18 

282.93 

-1.21 

.00 

299.79* 

-1.76 

.01 

19 

284.48 

-2.25 

.03 

301.71* 

-4.96 

-.16 

Ml 

219.98 

g** 

-.02 

253.41 

0** 

.04 

M2 

245.33 

-.26 

.01 

252.26 

.87 

.08 

M3 

226.44 

.19 

.01 

221.64 

1.62 

.01 

M4 

217.08 

.24 

-.06 

205.38 

1.95 

-.06 


*DAYTIME CASE 

**USED TO COMPUTE EMISSIVITY 



those found in Table 2.5 representing the quality of fit to the eight climato- 
logical profiles. 

The large differences in observed and ccraputed brightness ten^rtures for 
channels 11 and 12 in the mid-latitude profile, and channel 12 in the tropical 
profile, are most likely the result of errors in the extrapolation above the in- 
conplete water vapor profiles. Hie climatological O 3 profiles did reasonably well 
in accounting for observations in channel 9. The errors introduced in the 
brightness tenure ture by the rapid algorithm transmittance model are insig- 
nificant except for those channels v^ose absorption is due primarily to 
water vapor and ozone. These channels have not yet been introduced in the 
GLAS processing scheme. When they are introduced in the near future, the 
model for these channels may have to be nt^iroved. 
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Chapter 3 

Theory and ;^plicatic« of the Glas Physical Inversion Method 

3.1 Overview of the GIAS Physical Inversicxi Method 

■nie GIAS physical inversion approach to teirperature sounding from satel- 
lite observations is fundamentally different fron the statistical regression 
itethod used operationally at NOAA/NESS in that heavy reliance is placed on the 
ability to model accurately the instrumental response to atmospheric and surface 
conditions, while no use is made of statistical relationships between satel- 
lite observations and atmospheric tenperature profiles. The method involves 
starting with a guess set of atmospheric and surface conditions from which 
expected brightness teitperatures for the satellite observations are ccmputed. 
Then, iterative relaxation of atmospheric and surface conditions is performed 
according to the difference between observed and ccsnputed brightness tenpera- 
tures until sufficient agreement is reached. If sufficient agreement can not 
be obtained, no retrieval is produced for that location. 

The basic advantages of the method arise fron the ability to utilize 
first guess and other a-priori information to account directly for all factors 
affecting the observations, and to determine areas where retrievals should be 
regected. In addition, GIAS retrievals can be derived in regions of high 
topography, where operational retrievals \i^re not processed during PGGE. Most 
important of all, a ccnplete set of auxilliary meteorological parameters can 
be simultaneously determined v^ich are alljconpatible with each other and the 
observations. 

The GIAS processing scheme is comprised sequentially of the following 
elements: 
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1. Calibration of the data. 

2. Averaging of the data into fields of view. 

3. Attacihing initial guess conditions, T^(P), q^(P), Pg. 

4. Conputation of clear-colum radiances R®, from the initial guess 
conditions. 

5. Determination of the surface emissivity at 50.3 GHz, to 
coirpute microwave radiances. 

6. Determination of the iterative reconstructed clear-column 
infrared radiances frati potentially cloud-contaminated radiance 
observations, R^. When it is determine to be too clouc^ for 
reconstruction of clear column radiances, no retrieval is performed 
and steps 7-10 are skipped. 

7. Determination of T^, the iterative sea/land surface torperature. 

8. Calculation of I^, the clear-column radiances using the iterative 
sea/land surface toiperature and tetperature profile. 

9. Catparison of ^ and R^ for the tenperature sounding channels. 

10. If sufficient agreement between R^ and ^ is not found in st^ 9, 

calculation of t'^'*'^{P), the next (iterative) tenperature profile, and 
return to st^ 5 to ccnpute and continue the iterative proce- 

dure. Otherwise, the iterative procedure is terminated. 

11. Set quality flags. 

12. Calculation of clcud field parameters under all conditions. 

3.2 Preparation of the Satellite Data for Analysis 

The HIRS2 and MSU data used by GIAS is the level I-B data obtained frcm 
NOAA/EDIS. Dhcalibrated observations are given for each channel at each spot 
observed by the instrument, vhich is earth located and flagged. Calibration 
coefficients are provided with the data. The data is calibrated at GIAS using 

the coefficients with the appropriate algorithm described by Lauritson et al . 
C19793. 
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Figure (1) sinews the scan patterns of HIRS2 and MSU. In the processing sys- 
tem used in analysis of the data, one sounding, representative of an cirea of 
125 X 125 km, was perfonted at a spacing of roughly 250 x 250 km. This resolu- 
tion and spacing was chosen so as to be ccxtparable to that used for cperational 
terperature sounding by NQAA/NESS. Observaticxis at the 125 x 125 km resolution 
were generated by averaging of the HIRS2 data as shown in the figure. First 
the HIRS spots are blocked into roughly 250 x 250 km areas in groups of 6 x 
10, 6x8, 6x7, 6x5, or 6 x 3 depending cn the location in the scan array. 
These blocks are further breken into quadrants, as indicated in the 6 x 10 
array. In each quadrant, the spots are divided into two sets, one being the 
half ccxitaining the highest radiances for the 11 wn window channel H8 on HIRS2. 
These spots are indicated as vhite in the figure. Radiances for all spots in 
a given set are averaged together for each channel. Each 125 x 125 km area is 
new characterized by two sets of HIRS2 radiances for each channel. These i^o 
sets wilt he refevved to in the future as "two fields of view . " The differences 
in the radiances in these fields of view will be utilized in the cloud correo- 
tion algorithm to be described later, vhich estimates the effective clear csolumn 
radiance vhich would have been observed in the 125 x 125 km area if it were 
cloud free. For each field of view, an effective satellite zenith angle is 
defined as the angle vhose cosine is given by the average of the cosines of 
all the spot satellite zenith angles in the field of view. 

The 125 X 125 km ejaadrant in which the sounding will be performed is 
chosen as the one containing the waaunest field of view as measured by the HIRS 
11 window channel. The latitude and longitude of the centroid of this quad- 
rant IS taken to be the location of the sounding vhich vd.ll be performed. The 
HIRS2 ctoservations in this quadrant are also colocated vtfith ihe MSU observa- 
tions vhose centroid is closest to that of the quadrant, and vri.th initial esti- 


45 



SUB ORBITAL TRACK 





178 8 km 


Fxg. 1. Scan pattern for HIRS2 (small spots) and MSU (large spots). Clie retrieval is perfomed every 250 x 250 kra 
(solid line box) area in one of the 125 x 125 km (dashed line box) quadrants. HIRS spots in each quadrants are averaged 
into 2 fields of view according to the warmest and coldest radiances from the llym window channel. The closest MSU 
spot is assigned to the quadrant. 
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itiates of tertperature profile, humidity profile, and surface pressure obtained 
frcm interpolating a 6 hour forecast generated ty the GIAS GCM CKalnay-Rivas et . 
al., 1977? Falnay-Rivas and Hoitsraa, 1979] to the appropriate location and time. 

The preparation of data, v^ich r^resents steps 1-3 described above, 
associates in every 250 x 250 Ion area observed by the satellite a time, latitude 
and longitude, 2 sets of HIRS2 radiances and zenith angles, a set of MSU 
radiances and an MSU zenith angle, and model generated estinates of surface 
pressure, humidity profile, and teirperature profile for the selected qaadrant. 
This is sufficient information for analysis of the data. In addition, the 
area contains the HIRS2 radiances in the single coldest field of viotf in the 
250 X 250 hm area, as seen by the llym window, and the average HIRS2 radiances 
for all spots in the 250 x 250 hm area. These will be used in generating an 
effective cloud hei^t and cloud fraction for the entire area. 

3.3 Conputation and Significance of the Surface Emissivity at 50.3 GEiz 

Steps 4 in the processing scheme, dealing with the corputation of channel average 

radiances from atmospheric and surface conditions has been treated in detail in 

Chapter 1. The calculation of accurate radiances as a function of surface and 

atitDs^heric conditions is essential for all subsecjuent steps in the processing 

system. 

The emissivity of the surface in the microwave region is much more vari^le 
than in the infra-red. At 50 GHz, the emissivity is typically .45 - .65 for 
open ocean, increasing with decreasing tenperature and increasing wind speed, 
and .90 - .95 for land. Sea ice has an emissivity of .70 or more, depending on 
its history, and snow has an onissivity of .90 or less depending on the depth 
(Staelin, 1981). 

Given a tarperature-humidity profile and a surface tenperature and pres- 
sxare, the microwave brightness tenperature at a given zenith angle is much more 
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sensitive to the surface eraissivity than an equivalent infra-red channel would 
be [Kbmfield and SussTcind, 19773. For accurate calculation of briefness tettp- 
peratures in microwave channel 2 , vhich has about a 10% contribution frcm the 
surface, a precise knowledge of the micrcwave surface emissitivity is needed. 

This parameter in turn also provides inportant information about the surface 
properties. 

The micrcwave eiassivity is calculated frcm the 50.3 GHz channel, as part 
of the iterative scheme, according to 

e= R-i-/Tdr-R^iT-i(P^) 

[Tg - Rit3'Ci(Ps) (3.1) 

vhere Rj_ is the 50.3 GHz observed brightness tetiperature, Tg is the iterative 
surface teitperature, and T(P) is the iterative atmosperic teiperature profile 
used in the calculation of the vpward and downward microwave fluxes emitted by 
the atmosphere. Ihe transmittance functions are corrected for temperature, wa- 
ter vapor, and zenith angle as described in Chapter 2, but possible effects 
of liquid water attenuation on the 50.3 GHz channel are not accounted for. 

It is assimed that in areas containing significant atmospheric liquid water, 
the retrieval would be re^jected either on the basis of too nuch cloudiness or 
for non-oonvergence of the solution as explained later. Ihe emissivLty determined 
frcm the 50.3 QEIz channel observation is used, together with the iterative 
tenperature profile, and surfece terperature, to calculate brii^tness tenperatures 
for the other MSU channels. 

3.4 Accounting for Effects of Clcuds on the Infrared Observations 

The radiance observed in an otherwise homogeneous field of view, containing 
partial hcmogeneous clcud cover a, xs given, within a reasonable approximation 
ty 

% = <^x,CIO + (1 “ ®‘)Fx,dR (3‘2) 
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where Ri ^orn and Ri,cLR radiances vfiiidh would have been observed xf the 

fxeld of view were ccn$>letely clou<^ or clear, respectively. Coi^jutatxon of 
clear-Golumn radxances Ri^cLR can t>e done routinely as in ecjiation 1.1, hut 
ccRputation of Rx,CLD requires accurate Rncwledge of the optical as well as the 
meteorological properties of the cloud. It is more advantageous to be able to 
account for the effects of clouds indirectly than to have to model their 
radiative transfer properties. As shewn ty auith (1968) and Chahine (1974), 
an estimate of the dear-column radiance, R^, can be reconstructed from the 
observations in the two fields of view according to 

^1 = Ril ^t%,l - (3.3) 


Where is the observation for cdiannel i in the field of view j and n is 
given hy gi/(ot2"ai). Ihe fields of view are numbered in the sense that 
Rg,l > ^8,2' Chce n is detemined, dear-column radiances can be reconstructed 
from the observations by using equation 3.3. An estimate of h is obtained with 
each Iteration. 

It is seen from equation 3.3 that large values of n will tend to anpli^ 
noise in the observations and are, therefore, undesirable. In the other extreme, 
n = 0 inplies field of view 1 is clear and ri = -0.5 is taken vhen it appeairs 
both fields of viev cire clear. As shewn by Qiahine £1974] and Halem et al ., 
[19783 , n can be determined from the infra-red observationj as part of an 
iterative scheme according to 


nN = 


p N _ T3 
Ky 


(3.4) 


1^7,1 - K-7,2 

vhere Ry^ is the cenputed dear-column radiance for the 15pm surface channel, 
using the Nth iterative temperature profile. In this case, the scheme will con- 
verge provided only 4.3pm infrared channels are used for teamperature sounding 
in the lower troposphere. The rate of convergence increases with the difference 
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between the surface texperature and the cloud-tcp tenperatare. Under soras 
hi^ noise, low contrast ccxiditions, divergent solutions can occur in the 
sense that an overestinate of rF will cause an overestimate of the reconstruc- 
ted 4.3vnn. Clear-colunn radiances vhich, in turn, will yield an increased 
lower tropospheric tenperature, produce an increased value of and 

lead to an increased etc. 

When a lower trcpospheric-sounding ndcrcwave channel is available, such as 
channel M2, a superior method for determing n can be used, making the estimate 
of n less sensitive to guess errors and alleviating the need for use of 
15|iti channels, vhich are significantly affected by H 2 O and O 3 absolution, in 
cloud filtering. r\ is determined as in equation 3.4 but with the 4.3Un surface 
channel, 13, used instead of channel 7. The microwave channel is used to cor- 
rect errors in R]^ 3 ^ due to errors in the iterative tenperature profile. The 
error in determined frctn ecguaticai 3.4 is a result of either use of an incor- 
rec±. testperature profile to estimate the clear column radiance, ccxtputational 
uncertainties such as the effect of water vapor on the transmittance functions 
of channel 13, ofcservational eirrors in or errors in the assunption of 

only one degree of non-hcsnogeneity in the cxnibined fielcJs of view. The error 
in Rj 3 ^ due to a wrong tenperature profile can be well accounted for ty adjust- 
ing the cxmputed brightness tenperature (equivalent black bo(^ torperature) for 
channel M2 according to 


T 13 




(3.5) 


vhere '^M 2 and Tj^ 2 ^ the observed and calculated microwave briefness tenp— 
eratures, T^ 3 ^ is the calculated dear-column brightness tenperature for chan- 
nel 13, and T 13 is the corrected dear-column bri'^tness tenperature for chan- 
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nel 13. Iliis correction is based on the approximation that a bias on the itera- 
tive teinperature profile in the mid to lower troposphere will produce approxi- 
mately the same error in cotputed brightness temperature in infra-red and micro- 
wave channels sounding that portion of the atmosphere. Ihe corrected dear- 
column radiance for channel 13 is then given ty 

k'iS® = 813 ^ 13 “ + - ’’m2“3 (3.6) 

and n is now cxwcputed accxjrding to 

(R 13^-Ri3,i)/(Ri3,i-Ri3,2^ 

If the observations in the two fiel(3s of view are sufficiLently close/ most 
likely both fields of view are either clear or overcast. We discriminate these 
two cases ty coitparing T ^ 3 , the corrected clear coluim bri^tness teitperature 
for chcinnel 13/ to Tp 3 ^ 1 ^ the observed brightness tenperature for field of view 
1 . If T 8 “^ ^hd 1 ^^ > 4/ the fields of view are considered too 

clcxif^ to do a retrieval. In the other limit, if <0 and 1 “ 2 ^ 

< 1"C, nis taken as -0.5/ that is both fields of view are considered clear. 

Utilization of the micrcwave observation not only speeds up convergence 
tander all conditions/ tut st^ilizes the soluticai in the sense that a positive 
bias in the iterative tenperature profile in the lower troposphere will not, to 
a first approxiitation, cause an Increase in n The actual estimate of the 
clcud height and amount is performed after the final solution is obtained and 
will be described later. The micrcwave oorrecrticai was not enployed in chapter 
1, dealing with ccitparison of cxxputed and reconstructed radiances in cIcxk^ 
conditions, because the true temperature profile was known. 

3.5 Determination of Sea-Surface Teitperatures and Ground Tertperatures 

The mam factors inf luencmg the accxtracy of retrieved sea-surface 
temperatures from infra-red window observations are the effects of clouds and 
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humidity on the observations. Cjserational surface teniperature sounders, such 
as AVHSR, utilize very fine spatial resolution observations, of the order of 
1 Icn X 1 km, in an attenpt to firxi ccttpletely clear spots. The current analysis, 
using catibined infrared and micrcwave observations, does not require high 
spatial resolution nor the existence of clear spots for the determination of 
accurate sea or land surface tenperatures, vAiich are determined as an integral 
part of sounder processing systsn. The effects of clouds on the observations 
are accounted for by use of equations 4, 7 and 8, giving the clear column 
radiances for cJ.1 infra-red cJiannels. The two 3. 7jim window channels on 
HIRS2, whose brightness t^tperatures cire relatively insensitive to humidity, 
are used for determination of surface tatperature rather than the llym 
window, vhich has been used operationally on AVHRR. 

Given Ti^, the clear radiances, for the two 3. 7-jim window channels, 

18 and 19, are recx5nstrucrt:ed according to equation 4. At night, the reflected 
solar radiation term can be neglected from equation 1 and surface tenperatures 
are easily obtained from each channel as 




Ti(Ps) 1 

- (1 - Ei)RiNtT3_N(Pg)-jB^(TN)dT 

o 




The downward flux, R^^t, is approximated as 


(9) 


Ri% 




= 2 COS0 /B (T^)dT. 


(3.9) 


o 

This approximation is based on the assunption of an optically thin atmosphere 
and a Lambertian surface. In general, Tg^^S ^s,19 found to agree 
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with each other to 1“C, even under partially clou(^ conditxons. The surface 
teitperature, xs taken as l/ 2 [Tg^ 2 ,Q^ + 

During the day, the effects of solar radiation on the 3. Turn channels 
mist be accounted for in ctotaxning accurate surface tertperature retrievals 
from these channels. The solar radiation reflected off clouds in the fxeld 
of vxew has alreac^ been accounted for by the clear column radiance algorxthm. 
If additional clouds are xn the path of xncxdent solar radiance with cloud 
fraction c, the solar radxation strikxng the ground wxll be attenuated by (1-c) 
The solar radxation reflected off the clouds will not be seen by the xnstrument 
because of its narrow field of vi®^. The net effect is to reduce the solar 
radiation by a factor of (1-c). 

One can attenpt to account for reflected solar radxation dxrectly by sub- 
tracting (Pg)frcm and siabstxtuting the result into equa- 

tion 3.8. In the case of c=0, HiTj^''(Pg), the mean solar radxation across the 
channel traversing the path from the sun to the earth and back to the satellite 
can be well estimated as 2.1Gir x 10“^ B^C5600 K] cos0jjT^(Pg, 9gpp) where 
0H xs the solar zenith angle and the transmittance xs cxstputed at an effective 
zenith angle, 8 eff» vhose secant is given by the sum of the secants of the 
solar and the satellite zenith angles. The case of cj^ is eqpLvalent to an 
effective reflectivity p = p"* (1-c) , 

This procedure is xnpractical because of the uncertainty xn p^, even if 
c=0. If the surface xs Laitibertxan and the emissivity is known, p''^^, the 
directional reflectance, xs equal to (1 - ei)/Tt. Significant errors of up 
to a factor of 2 can be made in these estimations of p'±, vAiich may produce 
errors of up to 10®C xn retrieved surface teitperature. These errors arise 
fron uncertainties in and the non-Laitibertian character of the surface. 

The same uncertainties in e^, however, do not appreciably affect the 
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calculated thermal radiatxon. Esturated values of 0.85 over land and 0.95 
over ocean are used in the analysis. Eather than assume a value for p, Tg 
aM are solved for in an iterative manner, assuming only that p is the 
same for both 3.7jjn channels. 

For 3.7ym sounding channel x, vje can write 






N 


ATM,! 




=: BiCTg) + diH^- = 


(3.10) 


where atitospheric contribution to the calculated dear-column 

radiance, d^ = ~p / and is given ty = HxTi"(Ps)/Ti(Ps)* The 
left-hand side of equation 11, and consequently A^, is taiown in a given 
iteration. Assuming p^ and are the sane for both 3. 7 tin channels, 
one obtains the equation 


Bi(Tg)-aBj(Tg) = Ai - aA j = A (3.11) 

vhere a=Hj_VHj". This non-linear equation in one urikncwn, Tg, is solved 
iteratively according to 

-h’v/Tg^'''^ A 

e = 

-h'^TgM B^(Tg’^) - aB^(TgM) (3.12) 

e 

vhere v = (\^^ + Vj)/2. This procedure converges rapidly. Once Tg is de- 
termined, d is calculated frcm equaticjn (3.10). This provides a value of p vhich 
IS used in equation 1 to correct the 4.3-im channels for reflected solar 
radiation effects. 

The Iterative ground taiperature is used to ccnpute the estimated clear 
column radiances, for each channel. Tenperature profiling utilizes a 

ccaiparison of these radiances computed for the tenperature sounding channels 
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fron 1^ and T^(P) to the Nth reconstructed clear column radiances, The 

differences are used to update the tenperature profile. Over ocean, clima- 
tological sea-surface toiperatures have accuracies of better than 2®C and can 
serve reasonably well for the catputation of clear column radiances. Under 
conditions vflien it is felt that it is too clouc^ to retrieve sea-surface tatpera- 
tures of greater accuracy than climatology, the sea-surface tenperature is fixed 
at Its climatological value for the purpose of radiative transfer calculations. 
This decision is made only in the first iteration. Climatology is used if either 

“^IS “'^18,1 ^ 20”C, or [Tg - T nr.TM l > S^C, or both T 13 - Tig^i > 10®C and iTg - 
T rr.TM l >3®C, that is, either the reconstructed brightness tenperature is very 
far from the observed brightness tenperature, indicating a large cloud correction 
IS necessary, or the retrieved sea-surface tenperature differs significantly 
fron climatology, indicating a potential problem. If both indicators of a 
problem exist, the tolerance conditions made are more stringent. During the 
day, the same criteria are used with the llym window channel 8 . The sea- 
surface tenperature is fixed to climatology about 3% of the tme. 

The accurate a-priori hncwleged of the sea-surface toiperature is also 
used to indicate low level overcast vdiich may have been missed by the cloud 
algorithm test for overcast described previously. If Tg - T pt.tm < -3®C and 

/s/ 

1^18 - < 2 1 Tg - TqlmI, that is the retrieved sea-surface tenperature 

IS more than 3® colder than climatology and the difference is greater than half 
the difference in reconstructed and observed brightness tenperatxxres, the fields 
of view are considered to be full overcast with low level clouds. No retrieval 
IS performed under these conditions. The retrieval is also rejected if the 
final retrieved sea-surface tatperature differs frcm climatology ty more than 
5®. 
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3.6 5he Atmospherxc Teitperature Relaxatxon Equatxcn ” The Inverse Problem" 

Forms of the relaxation method of finding a solution to a set of radiative 
transfer equations was developed by Qiahine (1968, 1970) and Smith (1970). It 
differs from other methods m that it does not attempt, in any iteration, to 
find a best soluticxi to the set of equations (observatirais) but only to provide 
a set of parameters giving better agreement of observed and calculated radiances 
than obtained in the previous iteration. The iterative method is corputationally 
fast and stable. Moreover, to solve the inverse radiative transfer equation 
It is necessary first to put the equation in approximate linear form with 
coefficients vAiich are profile-dependent. Therefore, an "exact" solution must- 
be Iterative in any event. 

In the iterative relaxation scheme, modifications to the teaiperature pro- 
file are made according to the differences of observed and coiiputed radiances 
in temperature sounding channels weighted ty the sensitivity of observations in 
these channels to changes in tenperature in different parts of the atmosphere. 
Given and R^^, it is most convenient to coipare the reconstructed and com- 
puted brightness teitperatures, T^^ and T^^, because of an essentially linear 
relationship between changes in brightness tenperature with changes in tempera- 
ture profile. 

The differences between T^^ and T^^ are used both to cotpute the next 
iterative temperature profile and to decide vhen to terminate the iterations. 

The iterative procedure is terminated if the root mean square differences of 
T^^ and T^^ for the temperature sounding channels is not at least 5% lower than 
in the previous iteration. Also, a maximum of 10 iterations is performed. 

When the iterative procedure is ccstpleted, the retrieval is rejected as non- 
convergent if this RMS difference is not less than 1°C. In addition, the 
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retrieval is rejected as non-convergent if the conputed brightness ten^jerature 
for nacrowave channel 2, v^ich gives a measure of the mean troposcfceric tenpera- 
ture, differs frcxti the observation by more than l^C. *11115 channel/ which is 
not used directly in the solution of the relaxation equations, gives a good 
check on the consistency of the solution with an additional ctoservation. 

Micrcwave channel 2 is used indirectly through the cloud filtering algorithm as 
shewn in equation (3.6), Hie frequency of occurence of non-convergence is roughly 
independent of cloud cover, but a tendency occurs for rejection when the microwave 
observations with the largest zenith angle are used. 

In order to use the difference of reconstructed and conputed brightness 
tanperatures in the taiperature sounding channels to estimate the error in the 
Nth iterative tenperature profile, it is useful to look at the response of the 
brightness terrperature of a channel to changes in atmospheric tenperature pro- 
file. 


Prom equation (1.1), we see that to a good approximation, brightness tempera- 
tures for two closely related teirperature profiles, T(P) and T(P)+ 5 (P) , will 
differ by 


where 


Ti[T(P) + 6(P)1 -Tj^[T(P)] = /Wj^(P) 5(P) dln(P) 


dT, 


/ 


dR, 


Wi(P) = 


V 


dR, 


R-. 


V 


dT 


T(P) 



(3.13) 


f3.14) 


assuming all else remains constant. Figure (2) shows the weighting functions 
which, as defined in equation 3.14, relate the change in brightness teirperature 
to the change in atmospheric temperature profile, for the seven channels used 


57 



original 

OF POOR QUA^-^TY 



Fig. 2. Weighting ftmctions, for a U.S. standard atmsophere at nadir viewing, for the seven 
ECRS2-MSU channels used to deteimne the tenperature profile. These channels are used to estima 
lean temperatures for the layeis bounded by *. 
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in determining the atmosE^ieric temperature profile. 
It can be shown that 


JWi(P)dln(P) « 1 - Ti(Pg) (3.15) 

Therefore, to a good approximation, for two profiles differing by a constant, 

Ti[T(P) + 6] « Ti[T{P)l + 5 (1 - ■^i(Ps))^ (3.16) 

that IS, a constant shift of temperature profile throughout the atmosphere 
produces approximately the same change in brightness teirperature, reduced by a 
small amount if the channel sees the surface. Moreover, if we lock, for exanple, 
at channel M4 in Figure (2) , we see the brightness ten^erature in that channel is 
virtually independent of ternperatures above 10 mb and below 300 mb and is mostly 
dependent on tenperature between 40 mb and 130 mb. !Hierefore, we can attribute 
a difference, ^ between observed dear-column brightness temperature, 1 ^ 4 , and 
cotputed brightness temperature, ' 1 ^ 144 ' ^ ^ comparable difference in the true 
layer mean tenperature 130 mb and 40 mb, = /t(P) dlnP, and that of the Nth 
guess. 

In order to reduce the effects of noise and stablizes the solution, it is 
desirable to average the estimates of mean layer tenperature 

^■^N+l ^ ,-y^N _ (3^17) 


where is the mean value of Wi(P) in layer 1 . There is no need for a one to 
one relationship between the number of channels and number of layers. 

In the analysis, mean layer tenperatures for the ten layers shown in Figure 
(2) are used to determine the total temperature profile T^'*‘^{P). To insure 
uniqueness, a constraint is put on the system that 


qN+1 - -jG + I 
k=l 




N+1 


(3.18) 
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vAiere is a global mean ten^ierature profile, F]^(P) are enpirical orthogo- 

nal functions at 52 pressure levels extending from 1000 mb to 30 mb, given by 
the eigenvectors, with largest eignenvalues, of the covariance matrix of a set 
of glct»al radiosonde profiles, and are iterative coefficients. The co- 
efficients Aj^ uniquely define the solution. 

The L coefficients, A^"^^ can be determined for each iteration from the N+1 

k 

estimate of M mean layer temperatures, T^+1, and any other information we may 
have. If for example, we want to find the L coefficients that compute mean 
layer temperatures that agree to the M specified values best in a least square 
sense, then 

aN+1 = (f»f)-1f«AN+1 (3.19) 


where A is the vector of L coefficients, F is an M by L matrix with elements 
®'ik mean value of Fj^(P) in layer i, and A^N+1 ig the difference 
between and Tq^^^, the mean layer temperature of the global mean. 


A N+1 _ m N+1 _ m 
1 “ -^1 ^G,l 


(3.20) 


Because of uncertainties in the mean layer temperatures, increased stability is 
obtained by finding the coefficients vAiich minimize a combination of the diffe- 
rence between estimated and computed mean layer temperatures, on the one hand, 
and maximizes the likelihood of the solution on the other. The appropriate 
equation is given by 


aN+1 = [F'F + Oh]-1f»AH+1 (3.21) 

where H is a diagonal M x M matrix with being the inverse of the fraction 
of total variance arising from eigenvector i, and is a constant. In practice, 
equation 22 is used with six EOFs and ° = 5 x 10“^. 

It IS seen from Figure (2) that vhile very little detailed information is 
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contained about temperatures abcwe 50 mb, the observations in channels 2 and 
M4 are still quite sensitive to tenperature changes above this level. It was 
found that vAien the difference in the observed brightness for channel 2 and 
that computed frcan the first guess was more than 5°C, the retrieval was usually 
of poor quality at pressures less than 200 inb but satisfactory above 200 mb. 

Ihese retrievals were flagged as good only in the troposphere. Occurrence 
IS almost exclusive^.y over land at the high latitudes. 

3.7 Comparison of the Relaxation Method with Regression 

Ihe net result of equations, 3.17, 3.18, 3.20, and 3.21 is the iterative 
equation 

rpNtl = rjiN + ^ L ^-^1 ( 3 . 22 ) 

where is the 52 level global mean temperature profile, is the Nth guess 
temperature profile, (^-1^) is the difference between the reconstructed 
clear column brightness temperatures and those computed in the Nth iteration, 
is the matrix of weighting functions, defined by equation 15 in the Nth 
Iteration but normalized so that the sum of W over all channels equals 1 for 
any layer, and B is given by 

B = F(F'F + F'L (3~.23) 

where L is the matirx which produces layer average values from point values 
(e.g. T = LT, F = LF, W = LW) , and i is identity matrix. Ihe matrix BW is the 
fundamental interpolation matrix which produces a change in temperature profile 
given a difference between observed and corputed brightness temperatures. It 
is composed of two elements, the profile dependant weighting functions which 
contain the atmospheric^physics, and the statistical B matrix vAiich results 
from the constraints on the solution. Ihe term [B-I] arises’ from the 
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e}?pansion of the solution about the global mean and tenls to further stablize 
the solution under high noise conditions. This term would drcp out of equation 
(3.22) if in equation (3.18) were replaced by i.e., if the iterative teirpera- 
ture profiles were ccaistrained as an expansion about the Nth guess rather than 
about the glctoal mean. 

While the form of equation (3.22) is similar to that used in regression analy- 
sis, there are a number of significant differences. Foronost amoung these is 
that equation (3.22) contains the full physics of the problem, allows for the 
incorporation of initial guess information into the solution, for the iterative 
treatment of the effects of clouds on the radiances, and for the ability to 
determine vhether a solution can be found with satisfactory agreement of 
observed and cotputed radiances. 


3.8 Determination of Cloud Height and Mount 

Given a tenperature profile, an effective cloud height and percent cloud 
cover vSiidh match the outgoing long wave radiation can be determined by use of 
equations (1.1) cind (3.2). The cloud hexght is effective because of the possi- 
bility of the presence of multiple cloud layers in the field of view. The cloud 
fraction is effective both for the above reason, and because it is the product of 
the true cloud fraction and the cloud emissivity. It represents the cloud 
opacity in the field of view. If we assume a single cloud layer, whose emissivity 
is and vhose top is at with tenperature T(P<^), \^ere T(P) is the retrieved 

tenperature profile, then equation (1) can be evaluated at any assumed cloud top 
pressure to give 


\,cld^^c^ 


e 

ic 


CBiCT(P^)P.(P^) 


+ 


B (T)dT] 


(3.24) 


Using equation (3.2), an effective cloud fraction consistent with the assumed 
cloud pressure, \(Pc)/ be determined for channel i according to 
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(3.25) 


%,clr " % 

“i(Pc) - ! 

%,clr ■“ 

where % is the observation for channel i and Ri,clr calculated clear 

colunm radiance using the retrieved temperature profile. For any set of chan- 
nels, and u can be determined which minimizes the difference between the 
<±)served and ccnputed radiances for the channels. In the current analysis', two 
channels were used and « and P,-. were found such that **i(Pc) ~ ~ 

In the global retrieval program, one temperature retrieval is performed 
every 250 x 250 kra, and is localized in the 125 x 125 km quadrant having the 
field of view containing the warmest observed brightness temperature. Ohis 
quadrant is chosen because it is assumed to be the single least cloud contami- 
nated field of view. If cloud parameters are retrieved from the same quadrant 
that the temperature field is retrieved from, the general cloudiness would be 
systematically underestimated. To provide an estimate of cloudiness over the 
entire 250 x 250 km area, radiances frcm all four quadrants are averaged to be 
used in equation (3.25) for determination of cloud fraction, given a cloud height. 
The retrieved temperature profile from the clearest quadrant is assumed valid 
for the entire 250 x 250 km area and is used in equation (3.24) to estimate the 
cloud radiance as a function of cloud top pressure. 

To maxmize stability and minimize the effects of errors and uncertainties, 
the numerator and deimoninator of equation (3.25) should be maximized. Therefore, 
both the full overcast and the observed radiances should be as different as 
possible from the clear column radiance. For this reason, channels 6 and 7, 
the two 15Vm channels sounding closest to the surface, are utilized for cloud 
height determination, and the observations used to determine cloud height are 
taken from the single coldest of the eight fields of view in the 250 x 250 km 
area. The cloud height is allowed to be between 850 mb and the tropopause. 


63 



Using this cloud hex^t, the effective cloud fraction for Idle entire area is 
determined from the radiances of channels 6 and 7 averaged over the entire 
field of view. 

The cloud parameters obtained are effective in the sense that they repro- 
duce the outgoing longwave radiation but not necessarily the detailed cloudiness. 
Under multiple cloud layers, for exanple, a single intennediate cloud height 
would be found and the total cloud fraction woiiLd be underestinrated. 

If the retrieval performed in the quadrant with the warmest observation 
has been rejected, cloud parameters can still be determined in an identical 
fashion, but the initial guess is used in equation (3.25) to cotpute clear and 
clou(^ radiances rather than the solution. In the special case vhen-overcast- 
Icw level cloudiness has been found, 100% cloud cover is set at the pressure at 
\hich the guess tenperature is equal to the retrieved surface tenperature. If 
this pressure is greater than 600 mb, a second layer of clouds is locked for if 
the brightness taiperatures in the coldest field of viav are significantly 
lower than in the warmest field of view. The procedure is identical to that 
described previously, but it is assumed that there is ccsiplete overcast of the 
lower cloud decik throughout the entire 250 x 250 1cm field of view and the sur- 
face contribution to equation (1.1) is taiken to come from the lower cloud level 
rather than the earth. 

This procedure can provide cloud parameters under almost all conditions. 
However, approximately 20% of the time, it is epparent from the observations 
that partial cloudiness exists, and no consistent cloud height and cloud frac- 
tion can be determined from using channels 6 and 7 and equations (3.24) and 
(3.25). In such situations, a cloud field is not returned. We are currently 
investigating a procedure to determine the amount and height of clouds by 
minimizing the residual of conputed radiances for more than two tropospheric 
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sounding ciiannels CChahioie, 1982] to produce iitproved cloud parameters mder 
all conditions. 
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Chapter 4 

Results for January 1979 

4.1 Yield of Successful Retrxevals vs. Cloudiness 

440,000 retrievals were run for the period Jan 5 - Feb 2, 1971, with one 
retrieval in every grid area as shewn in Figure (1). Of these, 60% were acceptable 
retrxevals and 40% were rejected. The retrievals are brdeen dewn xnto quality 
flag vs. retrieved per-cent cloudxness in T^le 4.1, The column marked % xndxcates 
the fraction of all retrievals havxng the appropriate quality flag. Retrievals 
of each quality flag cure further broken down into the precentage of tiite they 
occur wxth a gxven cloud fractxon. The column marked * represents retrxevals 
in vAixdh no cloud hex^t and cloud fractxon could be retrxeved. This occured 
about 20.6% of the txme, most likely in areas contaxning more than one degree 
of homogeneity, either xn the sense of itultiple cloud formatxcns or varying 
scenes in the clear part of the field of vx®?. Retrievals of the first three 
types are all considered acceptable retrievals and statistics are also gxven 
for the total class of acceptable retrxevals. Retrxevals of the next four 
types were all rejected for the xndxcated reasons. Curaulatxve statxstxcs are 
also given for all rejected retrievals as well as for all retrxevals. The 
last row of the t^le, unlike the other rows, does not give the percentage of 
a given class of retrievals contaxnxng a given cloud amount, but xndxcates the 
percentage of all retrievals with a given cloud amount vhxch vere considered 
good retrievals. 

A total of 60% of all retrxevals were acceptdDle, 45% wxthout qualxficaticsi, 

3% with no sea surface toiperature retrxeved, and 12% with no stratospheric retrxe- 
val. For the class of accqptable retrxevals, 36.7% were found to be in clear 
cases, 21.7% were in cases vhere no cloud fxeld could be produced, 25.6% were 
xn cases with cloud fractions greater than zero but less than or equal to 40%, 

15% were xn cases of cloud fraction greater than 40% but less than or equal to 

PHECEDlf^G PAGE NOT FILMED 
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Table 4.1 

Per-Cent Retrieval Type vs. Per-Cent Cloudiness 
FRACTIOML CLOUD COVER 


Quality Flag 
good retrieval 


good retrie- 
val, no SST 
retrieved 


jcLR 1 .01-. 1 I .11-.2 1 .21-. 3 

I I ! i 

|40.8| 2.8 I 6.9 I 8.3 


31-.4 i .41-.5 ! .51-.6 i .61-.7 1 .71-.8 1 .81-.9 1.91-l.Ot 

I 1 I ^ ^ 1 1 

7.7 1 6.3 I 4.4 1 5.4 I .7 I .1 1 o I 


10.6 1 9.2 I 5.9 i 2.8 


1 — 
9. 

>6 1 

1 — 

1 3. 

1 

.4 

i 

2.8 1 
1 

1 

1 

1 

.8 ! 

1 




per cent good 
for given 
cloud fraction 


*no cloud field retrieved 
**acceptable retrieval 
***re]ected retrieval 




































































70%, and only 1% were in cases of cloud fraction greater than 70%. fbr the 
class of re 3 ected retrievals, 29.7% were clear, 18.7% had no cloud field returned, 
16.5% were xn cases of less than 46% clcudxness, 14.2% were in cases between 
40% and 70% cloudiness, and 20.9% were in case of cloudiness more than 70%. 
Success&al retrievals outnumber rejected retrievals for all cloud fractions up 
to 60%, and are of cotparable amount xn the 60-70% range. It is interesting 
to note cilso that for up to 70% cloud cover, the majority of rejected retrievals 
occur because of non-convergence, and that -Uie dxstributicn of retrxeved 
cloud amounts in non-convergent retrievals is sxmilar to that in accepted 
retrxevals. This xndxcates that successful HIRS2 retrxevals can be performed 
xn areas with cloud fraction of up to 70%. Retrievals indicating cloud fraction 
of over 70% in the 250 x 250 km area represent only 9% of the total cases. 

Figure (3) xndicates the distribution of successful retrievals xn terms of 
coverage per-day on the 4* lat. x 5'* Ion. grxd used xn the GIAS general circu- 
latxon model CKalnay-Rxvas et al . , 1977] . Each grid point is counted 1 or 0 
depending on whether a successful retrxeval is found in a gxven 6 hour perxod. 

The maximum possible yield is, in general, 2 retrxevals per day because most grxd 
points (except for high latitude points) eire covered only twice daily. As shewn 
xn Fig. 3, retrieval coverage of at least 1 time per day is almost global with 
the exception of a few areas vhere persistant cloudiness occured for moderate 
periods of time. Conspicuous by their ^sence are features due to the ITCZ or 
mountain eireas, indicating that these factors do not significantly decrease the 
retrieval yield. 

4.2 Tertperature Retrieval Accuracy 

Accuracies of the retrieved tearperature profiles for Jan 5 - Feb 2, 1979 
are shown in Figure (4), ccatparing mean layer tatperatures of GLAS retrievals, 
in the 9 pressure intervals shown, with mean layer teiperatures reported by 
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DAILY GLAS RETRIEVAL COVERAGE JANUARY 1979 
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Fig. 3. Average number of times per day at least one successful retrieval was performed per six hour period in a 
4° X 5° box. "Perfect" coverage is ncminally 2 times per day. 
















PRESSURE (MB) 


RMS LAYER MEAN TEMPERATURE ERRORS COMPARED 
TO RADIOSONDES ±3 HR ±110 KM 
JANUARY 5 - FEBRUARY 2, 1979 


RMS ERRORS BIAS ERRORS 

GLAS RET (8150) 2.20 QLAS 

NESS RET (5486) 2.41 NESS 



DEGREES (K) DEGREES (K) 

Fig. 4. tarors of retrieved Dean layer temperatures compared to colocated radiosonde for 
January 5-February 2, 1979. 
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radiosondes colocated in space to 110 l<m and in tine to 3 lir. Retrievals flag- 
ged lead are not xncluded in the statisitics. In addition, retrievals flagged 
bad in the stratosphere are not included at 200 itib and above. Also shown are 
retrieval accuraacies of the operational NESS retrievals fdor the same period 
of tine (subsequent to this period of tone, NESS introduced changes to their 
cperational clear column radiance algorithm CMcMillin and Dean, 1982], but 
this new processing syston has not been applied to the Jan. 1979 data) . In 
the case of NESS retrievals, reported mean layer virtual tertperatures are 
cctipared to virtual tenperatures derived from the radiosonde tenperature- 
humdity profiles. The GLAS retrievals cire seen to be significantly more 
accurate in the troposphere, thou^ sli^tly less accurate in the stratosphere. 
The total RMS error of the 8150 colocated GLAS retrievals is 2.20®, .21® lower 
than that of the 5486 colocated NESS retrievals. Also shown in Figure (4) are 
the bias errors of the retrievals. The GLAS retrievals tend to have a .4® oold 
bias in the troposphere and a .4® warm bias in -the stratosphere. Ihe causes of 
these biases are being investigated. A tendency has been found for the bias to 
disappear in clear areas and be largest in tropical clou^ areas. 

Figure (5) shows a breakdown of the error statistics into categories of 
varying cloud cover for the period Jan 5 - 15, 1979. Ihe NESS retrievals were 
stratified according to the reported retrieval type. In the retrievals marked 
clear NESS treated the HIRS2 observations as not cloud contaminated and they 
applied no cloud correction to the HIRS2 radiances. In the retrievals marked 
N*, NESS performed a correction to account for cloud effects on the HIRS2 obser- 
vations before the retrieval was performed. In the retrievals marked clouc^, 
the effects of clouds on the HIRS2 observations were considered by NESS to be 
too large to be accurately accounted for and only HIRS2 channels 1-3, sounding 
the stratosphere, were used in the retrieval together with the MSU and SSU 
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OF POOR QUALITY 

MEAN LAYER TEMPERATURE ERRORS COMPARED 
TO RADIOSONDES ±3 HR +110 KM 
JANUARY 5-15, 1979 


GLAS (2647) NESS (1646) 


a < .1 (2.21) 57% CLEAR (2.25) 71% 

.1<a<.4 (2.18) 23% (2.44) 9% 

« > .4 (2.35) 20% CLOUDY (2.86) 20% 



2.0 2.5 3.0 2.0 2.5 3.0 3.5 

DEGREES (K) DEGREES (K) 

Fig, 5, BreaJsdown of RMS mean layer tonperature errors vs. degree of cloudiness for January 
5-Janijary 15, 1979, GLAS retrievals are separated according to retrieved cloud fraction, a. NESS 
retrievals are separated according to reported retrieval type. 
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observations. Hie GIAS retrievals were partitioned according to almost clear 
a < .1)/ partially clou<^ (.1 _^a ^ .4), and lii^ly cloucty (a > .4) conditions. 
The latter cutoff appears to be about the region vftiere NESS began to apply 
their cloudy algorithm. NESS has a considerably higher percentage of "clear" 
retrievals than GLAS has. The main reason for this difference is that the 
GLAS cloud fraction refers to the entire 250 x 250 km area vSiile NESS clear 
refers to the existence of clear 30 km spots. 

The accuracy of the GLAS retrievals is seen to degrade much less with 
increasing cloudiness than that of the NESS retrievals. GLAS retrievals with 
cloud fractions between .1 and .4 are in fact quite cotparable in accuracy to 
those obtained under clearer conditions. The NESS N* retrievals shews a large 
degradation over the clear retrievals in the upper troposphere, possible due to 
the effects of multiple cloud layers on the radiances. Such cases would hope- 
fully be identified and flagged in the GLAS retrieval system. The NESS clear 
retrievals appear to degrade scxnevhat in the 400-500 itib layer, possibly also 
due to a residual cloud effect. The NESS clou(^ retrievals, using only 2 micro- 
wave channels to sound the treposphere, are of significantly lower quality, pri- 
marily due to lack of data. The GLAS retrieval system shews that reasonably 
accurate retrievals, using both the HIRS and MSU channels, can be performed 
under almost all cloud conditions. 

Figure (6) shesra the RMS errors in the thickness between the mandatory 
levels and 1000 nib obtained by cewparing thicknesses ooitputed from the retrievals 
and those reported by the colocated radiosondes fqr the same period as m 
Figure (5). The ratio of a thickness error to the thickness is rou^ly propor- 
tional to the error in the average retrieved tenperature in the layer to the 
average temperature in the layer. If the average retrieved tenperature error 

J 

were constant, we would expect the thickness errors to grew linearly with the 
log of the pressure. 
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GLAS RETRIEVAL RMS THICKNESS ERRORS 
COMPARED TO RADIOSONDES i3 HR ±100 KM 
JANUARY 5-15, 1979 



/IZ METERS 

Fig. 6. RMS error of retrieved thickness between pressure P and 1000 nib coipared to colocated 

reports, broken up according to retrieved cloud fraction, a. January S^anuary 
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The thickness errors in the GLAS retrievals appear to grow less than 
linearly in the log of the pressure. This inplies the error in the average 
tenperature between the surface and hei^t z is decreasing with hei^t, vhich 
is reasonable because of possible cancellation of errors in the detailed tstpera- 
ture profile. The RMS errors do not degrade appreciably Vidth increasing cloud 
fraction, ranging frcra 33-43 meters at 500 mb, 51-59 meters at 300 itb and 
65-84 meters at 100 nto. In terms of average thidkness taiperature errors, 
these values correspond to 1.6-2.1®C at 500 iib, 1.3-1.6‘’C at 300 mb and 1.0-1.2°C 
at 100 itib. The average errors- decrease with height and-are all considerably 
less than the mean layer errors shewn in Figure (5). 

4.3 Sea Surface Teitperature Field for January 1979 

The sea/land surface tenperatures produced the GIAS retrievals can be 
used to produce global rronthly mean fields of tenperatures and their diurnal 
variations. In particular, monthly mean sea-surface toperature fields are very 
significant for climatological studies. Conventional in-situ sea-surface ton- 
perature measurements frem ships and bucys are numerous in the northern hemi- 
sphere but coverage is sparce in space and time in the trcpics and southern 
hemisphere. Sea or land surface tenperatures, averaged over the 125 x 125 km 
area, are retrieved fron each successful HIRS2/MSU retrieval except those in 
vhidh the sea surface teiperature was held fixed at climatology. At ni^t, 
equation 3.8 is used to obtain the surface tanperature v^le eqpaations 10-12 
are used during the day. Figure (7) shows the mean sea-surface teiperature 
field derived for the period January 5 - February 3, 1979, obtained ty avera- 
ging all sea-surface teiperature retrievals in 4® ty 5® latitude-longitude 
bins, tfo smoothing, rejection criteria other than those described in the 
text, or adjustments for bias removal, were applied to the data. Differences 
between SST analyses using only ni^t (3AM local time) retrievals and using 
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MONTHLY MEAN SST JAN 1979 FROM HIRS2/MSU 








20W SOW 40W 0 40E 80E 120E 160E 

Fig. 7. Monthly mean sea-surface temperature field for January 1979, obtained by averaging uH successful sea surface 
temperature retrievals in the 4® Lat x 5° Lon grid boxes and contouring the results. No empirical correction or filterir^ 
ms done to the data. The surface tarperature is determined primarily by the 3.7ym and 4.0jjm window channels. 
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only day (3PM local tame) retrievals were very small and never mere than 2® in 
the open ocean. This indicates that the procedures used to account for the 
effects of solar radiation on the 3.7 tm and 4.0 pa diannel olDservations can 
produce accurate daytime sea-surface tarperatures using these short wave 
window channels. Both day and ni^t retrievals are included in Figure (7). 

The major climatological sea-surface terrperature features, such as the Gulf 
Stream, the I^aroshio, Humboldt, and Benguela currents, and the sea surface 
tenperature ininimum at the equator in the eastern Pacific are readily observable. 

Of particular interest is the sea-surface tenperature anomaly field for 
January 1979 shewn in Figure (8a) obtained by subtracting the NCAR climatology, 
based on an average of data from 20 Januaries, from the GIAS January 1979 sea- 
surface tenperature field. The deviations from climatology are small, being 
less than 2® in all cases with the exception of the centers of the warm anomaly 
west of South America, and the cold anomaly midway between South America and 
Australia. 

Figure (Sb) shows the ancmoly field for January 1979 computed by subtracting 
an analysis based on ship and buoy measurements, coitpiled by Fleet Numerical 
Weather Center, from the same NCAR January climatology. The Fleet Nimerical Anal- 
ysis can be taken as a measure of ground truth in the areas of dense coverage. 
Agreement of the major anerraly features in the Northern Hemisphere, such as 
the cold Pacific areas, north of 40®, off the west coast of North America, and 
centered at 180°W, 10®N, as well as the warm Atlantic off the west coast of 
Africa, is excellent. In looikmg at this map, it should be remstibered that no 
bias errors were ranoved from the retrieved sea-surface torperatures . Even 
small biases of a few tenths of a degree would have a significant effect on 
the location of the contour lines, especially the 0® bias line, vhidi also 
matches extremely well. We can conclude that the absolute accuracy of the 
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Fig. 8a. Difference between the sea surface temperature field in Fig. 7 and a climatology based 
on an average of 20 Januaries conpiled by NCAR 
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Fig. 8b. Difference between the sea-surface taiperature field for January 1979, canpiled by Fleet 
Numerical Weather Central fron ship and buoy roeasursnents, and the saute NCAR January climatology. 


SHIPS AND BUOYS SST 




79 




climatolgical sea-surface tenperature data is quite high. Detailed coiparison 
of the analyses gave RMS differences of .4®C in the North Atlantic Ocean and 
.6® in the North Pacific. 

No significant difference is found in the two fields in the trcpics. 

This contrasts to the findings of Barnett et al . / [1979!], with regard to opera- 
tional sea-surface teaiperatures in the central trqpical Pacific Ocean, which 
are shown to have very large biases (l®-4°) cxstpared to XBT raeasuronents/ with 
errors that are highly correlated with cloudiness and huiradity. They concluded 
that the spatial structure of the error field, and the large magnitude of the 
errors, preclude the use of cperational sea-surface tenperature for climato- 
logical studies of the area. Mc:Clain (1981) has shewn that inproved sea-surface 
terrperatures can be retrieved fran the more recent 5 channel version of the 
AVEIRR flying on NOAA 7 and the cperational procedures have new been modified. 

Agreement in the Southern Hemisphere is also extremely good. Note for 
exanple the excellent agreement of the oscillating warm, cold, warm ancamly 
pattern in the latitude band frem 10®S to 30°S. The largest difference in the 
fields occurs south of 40“S, vhere the conventional data field is noisy as a 
result of sparce data. The apparent large difference in the fields between 
60®E-120®E and 40®S-50°S is in fact only a 1® difference in teitperatures in 
this area. In general the Southern Hemisphere ancmaly field retreived from 
the satellite data is less noisy and better defined than that frem conventional 
data. 

The eastern trepical Pacific, between 10®N and 10®S, has basically no 
anomaly in either field. An extended warm ancamly in this area is indicative of 
the "El Nino" phenonenon, vhich has been associated with major effects on the 
Climate of the Northern Hemisphere CRasitusson and Carpenter, 1982; Horel and 
Vfellace, 1981]. 
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Detailed analysis of specific sea and land surface teitperature fields and 
coitpariscns witli colocated surface measxnrenents mil be treated in a separate 
publication. 

4.4 Effective Cloud Cover and Cloud Hei^t for January 1979 

Effective cloud heists and cloud fracticxis, consistent with up^lling 
longwave radiation, are derived for each 250 x 250 km area through use of equa- 
tions (3.24) and (3.25). In the special case of lew level overcast, 100% cloudi- 
ness is assigned to the level v^ose tesrperature equals the retrieved surface 
'tenperature. The cloud top pressure and cloud fraction are effective for two 
reasons. First, possible multiple level cloud formations are assigned a single 
effective cloud top pressure, which should lie scsnevfliere in the range of the 
cloud top pressures in the field of view, generally close to that of the highest 
(coldest) clouds. Second, the effective cloud fraction is a functicn of the 
conputed effective cloud top pressure for a given situation, with decreasing 
cloud tep pressure (increasing cloud top height) corresponding to descre^ing 
cloud fraction. Consequently, cloud cover frem lew level clouds may be signi- 
ficantly underestinated. Moreover, even for single layer clouds, the effective 
cloud fraction corresponds to the cloud onissivity times the actml cloud 
fraction. This means that cirrus clouds mil be underestimated ty up to a 
factor of 7 and mid level clouds ty up to 30% CPaltridge and Platt, 1976; 

Chahine, 19823 • 

Monthly mean fields of cloud fractions, cloud top pressures, and cloud 
top tsnperatures were produced ty averaging the retrieved cloud parameters 
in the 4“ x 5“ grid for the period of Janaury 5-February 2, 1979. The 20% of 
the cases vhere no cloud field was retrieved was not included in the averages. 
The average cloud top pressures and tenperatures were taken as the average of 
appropriate quantities wei^ted ty the corresponding cloud fractions. Cases 
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in vSiich no cloud field can be retrieved tend to have Lcw-intennediate total 
cloud cover and their omission most likely doe not significgintly effect the 
average cloud statistics. 

Figures (9a-c) show ccxxtours of average cloud top fractions, cloud 
top pressure, and cloud top teirperature for the period. Ihe major features 
such as the Intertrcpical Convergence Zone, the storm tracks over the Atlantic 
and Pacific Oceans, the oceanic deserts, and the region of the Siberian high, 
are clearly visible in the cloud fraction itap. In the cloud top pressure map, 
in vfliich only areas of high {<300 mb) and lew (>600 mb) clouds have been indi- 
cated, the Intertrcpical Convergence Zone is again clearly visible as well as 
the extensive areas of predominantly stratus cloud cover below 50 "S and off 
the west coast of ^ the southern hemisphere continents south of 20 “S. Ihe global 
average effective clcud cover obtained for January 1979, sanpled at 3 A.M. and 3 
P.M. local tune, is 25%, a value considerably lower them the accepted value of 
50% global cloud cover based on ground observations [Paltridge and Platt, 1976]. 

The clcud tep tenperature fiel d appears similar to the cloud pressure map 
in the tropics but very cold mid-level clouds are found in the northern hemi- 
sphere over land. It is interesting to note that the -40 "C contours, corre- 
sponding to the formation of ice clouds, closely follows the continental con- 
tours in the northern hanisphere. 

Verification of retrieved cloud parameters is difficult. Another repre- 
sentation of cloudiness for the same period is shown in Figure (10), [Gruber 
and Vamadore, 1981], containing the averaged outgoing longwave radiation, as 
determined from the Hum channel on AVHRR on TIROS-N, for the period January 7 - 
January 31, 1979. In the trcpics, outgoing longwave radiation is a direct cjuali- 
tative measure of cloudiness because, as shewn in Figure (9), areas of extensive 
cloudiness contain cold clouds, vhile clearer areas contedn warm clouds or allow 
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GLAS CLOUD FRACTIONS JANUARY 1979 



Fxg. 9a. Average effective cloud fraction retrieved by GLAS for January 1979 using coribined 3 A.M. and 3 P.M. soun 
dings . Ihe cloud parameters are detennined primarily from channels 6 and 7 , the two 15ym channels sounding lowest 
in the atmosphere. 
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RETRIEVED CLOUD TOP TEMPERATURES JAN 1979 
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Fig. 9c. Average effective cloud top tonperaturej weighted by cloud fraction, retrieved by GLAS for January 1979 
Tenperatures are in “C. a 








radiation from the vrartn surface. Tlie najor patterns of cloud cover inferred 
frcm AVHRR/ and detennined from the HIFS/MSU system, agree almost perfectly in 
the trcpics. In the e3Ctra*-tr epics, especially in the northern heinisjiiere winter 
over land v^ere cold land tenperatures exist, and, as shewn in Figure (9), 
cold, thou^ not necessarily extensive, clouds exist, outgoing longwave radia- 
tion is not necessarily a measure of cloudiness. Note for example, from 
Figure (9) that the North Atlantic Ocean is considerably more cloudy than either 
North America or Eurasia. The outgoing longwave radiation in this area, is 
higher than that over the continents, however, as shewn in Figure (10). 

4.5 Micrciwave Surface Emissivity - Ice and Snew Cover for January 1979 

The emissivity of the surface in the microwave region is a strong function 
of surface conditions. At 50 GHz, open ocean has emassivity values ranging from 
.45-. 65. The emissivity increases with decreasing teiperature and also increases 
with increasing foam cover vhich is a measure of wind speed CWilheit, 1979]. 

Land has amissivities typically greater than .9. Ice over ocean would have an 
emissivity of the order of .7 for old ice, .9 for new ice, and .7 or more for 
mixture of n«v ice and open ocean. Snow over land would have an omissivity of 
the order of .90 or less [Staelin, 1981]. Thus, the surface emissivity can 
give a measure of snew cover land, ice cover over ocean, and possibly also 
boundary layer windspeed over epen ocean. 

Passive microwave sounders designed to measure surface properties, such as 
ESMR and SMMR, have used frequencies less than 40 GHz so as to avoid attenuation 
by atmospheric oxygen, vhich mixes atmospheric effects into the signal. Use of 
equation (3.1) accounts for atmospheric effects, and also allows for the incor- 
poration of surface temperatures, obtained from the infra-red 3.7pm channels, 
into the determonation of surface emissivity from the HIRS2/f®J system. The 
spatial resolution of MSU is considerably lower than that of ESMR or SMMR 
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JANUARY 1979 OUTGOING LONGWAVE RADIATION (WM'^) 
FROM AVHRR (DERIVED BY NOAA/NESS) 
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Fig. 10. Oitgoing longwave radiation field for Janmry 1979 derived by NOAA/NESS from the llvm window AVHRR data 
on 3 A.M. and 3 P.M. observations. 







"however/ and this will degrade hi^ resolution features vhich could be sxgnifi- 
cant for some applications. 

Figure (11) shows the surface emissivity/ averaged over the 4® x 5® grid 
for January 1979. As esipected/ the most o"bvxous features eire the continents/ 
showing rapid gradients from .9 - .6. The intermediate entassivity values are 
partiadly due to ccxitouring, and partxally due to soundings with mixed fields 
of view at the coastlines. The emissivity'in South America/ Africa, and Australia 
is uniformly greater than .9. North America/ Eurasia/ and Antarctica shew 
ccoplicated patterns containing lower emissivities, indicative of snow cover. 

Ihe oceans show the generally espected features of emssivity "between .45 
and .65/ with eraissivity increasing with decreasing sea surface tenperature. A 
distinct .7 contour, indicated "by the solid line in the Figure (11), is observed 
in the Antarctic Ocean and in the Bering, Labrador, Greenland, and Barents 
Seas. We interpret this contour to be a measure of the sea ice extent in 
these areas. We have no verification of sea ice in the Antarctic "but do have 
a sea ice field produced fran SMMR in the northern hanisphere CCavalieri, 

1981]. Figure (12a) shews the sea ice line as deduced from HIRS2/MSU by the 
.7 enissi"vity contour over water. Figure (12b) shews the sea ice line determined 
from SMMR, with a 25 km resolution. The ice margin from the ®1MR data was 
taken as the contour of sea ice coverage greater than 40%. The agreement is 
quite good, considering the difference in resolution of the instrumente. 

The details of the emissivity over land may be indicative of snow cover. 
Figure (13a) shews the averaged surface anissivily over Nosrth America. The .7 
emissivity contours^ indicative of the onset of significant amounts of sea ice, 
and the .9 contour, indicative of typical snow free land, are marked. Most of 
North America has land emissivities less than .9. Figure (13b) shows a itap of 
monthly mean ice and snow cover derived "by averaging the weekly observations 
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SURFACE EMISSIVITY JANUARY 1979 
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Fig. 11. The averse global surface smssivity at 50.3 GHz derived by GLAS for January 1979. Resolution is at the 
4® X 5® grid. Ehussivities for all the zenith angles of MSU have been averaged together without correction for angle. 
The .7 emissivity contour over ocean is indicated by the dark line. This contour is taken to be indicative of sea-ice 
extent . 




HIRS2/MSU ICE EXTENT JAN 1979(125 KM) 
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Fig, 12a. North Polar projection showing average sea-ice 
extent for January inferred frcxn the .7 emissivity contour 
shown in Fig. 11. Sea ice is indicated by the dots; open 
ocean by the waves. The MSU spot size is 125 x 125 km 
at nadir. 


Fig. 12b. Sea ice extent for January 1979 inferred by 
CJavelieri (1981) from analysis of SMR data, with a 25 x 
25 km resolution. Tbe ice/ocean line was taken as the 
line which represented 40% ice cover for the period 
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SURFACE EMISSIVITY JANUARY 1979 



Fxg. 13a. Blowup of the northwest quadrant of Fig. 11. Ihe .7 anissivity contour over ocean, indicative of sea ice 
extent, is shown as the solid line, as in Fig. 11. The .9 contour, indicative of land is shown as the dashed line. 
Near ocean, it represents the land boundary. Inland, as in North America, it is indicative of the onset of snow~ 
covered land. Further north, it is indicative of solid newly frozen sea ice. 
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PERCENT SNOW COVER JANUARY 1979 


Fig. 13b. Itercent snow cover for January 1979 obtained from the NOAA data of Dewey and Heua (1981). The contours indicate 
the aresis of 30%, 60% and 90% snow cover for January 1979. 
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produced cperationally ty NQAA/NESS CDewey and Heim, 19811 xnto the 4® x 5® 
grid toxes for the nonth of January, The .9 emissivity ccsitour lies oonpletely 
in the snow firee area while the .85 contour closely ^proxuitates the 60% sncw 
cover contour, indicated in Figure (13b). Detailed studies will be conducted 
to quantiJ^ the relationship between snow cover and ice cover and surface 
emissivity at 50.3 GHz, including the ability to distinguish between old and 
new ice and to estimate snow depth. The ability to get qualitative sncw and 
ice sncw maps frcsti the operational TIROS-N sounding system is well demonstrated 
by the current results, 

4.6 Summary 

The GIAS processing system for analysis of HH^2/MSU data by finding 
atmospheric and surface conditions vhich are solutions to the multispectral 
radiative transfer equations has been shewn to produce not only atmospheric 
temperature profiles of hi^er accuracy than those produced cperationally during 
the same period, but eilso to produce monthly mean fields of surface torperature, 
cloud prcperties, and ice and sncw cover vhich shew reasonable agreement with 
ground truth. Simultaneous analysis of infra-red and microwave observations 
erhances the strength of observations in each spectral region. Infra-red 
soundings have higher sensitivity to Icwer trcposphenc testperatures and ground 
teirperatures. Micrewave obseiivations are better for sounding the stratosphere, 
are sensitive tO snow and ice cover, and greatly aid in correcting the infra-red 
observations for cloud effects. 

The retrieval system has the potential of iitproving the iirpact of HIRS2/MSU 
data on mid-range forecasting (Halem et al . , 1982) both because of inproved 
accuracy of the retrievals and also because the retrieval process can be coupled 
directly into the forecast/assimilation cycle, using the current forecast as a 
first guess. This provides an analyzed field vhich is consistent not only 
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with cxxiventional observations tut also with satellite radiance observations 
as well. The itonthly mean fields can form the basis for a long term olirnatology 
data base, derived froti sinnilar instrumentation, starting from 1979 and continuing 
throu^ the 80' s. Modifications to the processing system are currently being 
made, based on subsequent research and ej^rience gained frcm cinalysis of the 
products of the current system. The inproved syston, as well as new results, 
will be described in a future publication. 
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